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Abstract 



Black holes have their proper oscillations, which are called the quasi-normal modes. The proper 

oscillations of astrophysical black holes can be observed in the nearest future with the help of 
gravitational wave detectors. Quasi-normal modes are also very important in the context of testing of 
the stability of black objects, the anti-de Sitter/ Conformal Field Theory (AdS/CFT) correspondence 
and in higher dimensional theories, such as the brane-world scenarios and string theory. 

This dissertation reviews a number of works, which provide a thorough study of the quasi-normal 
spectrum of a wide class of black holes in four and higher dimensions for fields of various spin 
and gravitational perturbations. Wc have studied numerically the dependance of the quasi-normal 
modes on a number of factors, such as the presence of the cosmological constant, the Gauss-Bonnet 
parameter or the aether in the space-time, the dependance of the spectrum on parameters of the 
black hole and fields under consideration. By the analysis of the quasi-normal spectrum, we have 
studied the stability of higher dimensional Reissner-Nordstrom-de Sitter black holes, Kaluza-Klein 
black holes with squashed horizons, Gauss-Bonnet black holes and black strings. Special attention is 
paid to the evolution of massive fields in the background of various black holes. We have considered 
their quasi-normal ringing and the late-time tails. 

In addition, we present two new numerical techniques: a generalisation of the NoUert improve- 
ment of the Frobenius method for higher dimensional problems and a qualitatively new method, 
which allows to calculate quasi-normal frequencies for black holes, which metrics are not known 
analytically. Also we considered a possibility of construction of the acoustic analogue of the Schwar- 
zschild black hole. 
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Chapter 1 
Introduction 



The general relativity implies two qualitatively new phenomena, that cannot be described within 
Newton's gravity. First, there exist objects with masses so high, that the escape velocity from them 
exceeds the speed of light. These objects were called black holes. Second, the general relativity 
predicts gravitational waves, which appear due to finiteness of the gravitational interaction speed. 
These waves may be detected in the nearest future with the help of gravitational antennas. During 
last years there have been an advance in construction of the gravitational wave detectors, such as 
LIGO, GEO, VIRGO, TAMA and the planned space-based detector LISA, and the detection of 
gravitational waves is expected very soon [1]. 

One of the most promising sources of gravitational waves is collisions of black holes and/or of a 
star and a black hole. The result of these processes is a black hole with higher mass, which absorbs 
the gravitational waves. Therefore, the gravitational waves quickly decay and, at sufficiently late 
time, can be considered as small perturbations of the black hole metric. It means that one can 
study the linear perturbations, neglecting the higher order corrections. This approximation for the 
late-time behavior of the gravitational waves provides a good accuracy, which was checked by full 
non-linear simulations of collisions of two black holes 

It turns out that the late-time behavior of the gravitational perturbations does not depend on 
the way they were induced. At the late-time stage of the evolution of gravitational perturbations we 
observe the damping oscillations, which give way to asymptotic tails at very late time. The damping 
oscillations are characterised by complex frequencies, which are called quasi-normal modes. The real 
part of a complex frequency describes the actual oscillation frequency, while the imaginary part is the 
decay rate of the particular oscillation. The set of the quasi-normal frequencies form the spectrum. 
Being dependent only on the black hole parameters, the quasi-normal spectrum appears to be an 
important characteristic of a black hole, or as it is said its footprint. Therefore, detection of the 
quasi-normal modes allows us to determine the black hole parameters and compare them with those 
obtained by astronomical expectations pj. 

Despite they are well described by linear approximation, the quasi-normal modes, when detected, 
can be used to check the general relativity as a full non-linear theory. This allows to study some 
aspects of the gravitational theory that cannot be experimentally confirmed without considering 
the regime of strong gravity. A good example is a phenomenon of violation of the local Lorentz 
symmetry. If we assume that there exist the locally preferred state of rest at each point of the 
space-time, we can describe this state by a unit time-like vector field. This vector field is called 
"aether" |1] . The parameters of the theory of the aether can be bounded by post-Newton corrections 
found from astronomical observations. Nevertheless, these bounds do not answer the question if 
the aether exists or not. Thus, one of the possible ways to observe the effect of the aether is the 
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determination of the corresponding shift in the quasi-normal spectrum of black holes. 

Another motivation for study of the black hole perturbations is checking of stability of the black 
holes. This is a very important property for higher dimensional theories, such as the brane-world 
scenarios and string theory [5] . Since in higher than four dimensions there is no uniqueness theorem, 
stability may be the criteria which will select physical solutions among a variety of "black objects": 
black holes, black branes, black rings etc. It is easy to understand that the quasi-normal spectrum 
of a stable configuration contains only damping modes, while black objects, unstable under small 
perturbations, must contain at least one growing mode in their spectrum. 

In addition, the quasi-normal spectrum can be interpreted in the context of the anti-de Sit- 
ter/Conformal Field Theory (AdS/CFT) correspondence |6]. A black hole in the anti-de Sitter 
space-time corresponds to the thermal state in the dual Conformal Field Theory. The tempera- 
ture of the thermal state coincides with the Hawking temperature of the black hole. Within the 
AdS/CFT duality the quasi-normal spectrum of the black hole in AdS corresponds to the poles 
of the retarded Green functions on the CFT side [7]. Due to the AdS/CFT correspondence we 
are able to calculate non-perturbative effects of the finite-temperature field theory at strong cou- 
pling by studying black holes in the anti-de Sitter space-times. The hydrodynamic parameters of 
the quark-gluon plasma were calculated within this approach, showing a good agreement with the 
results obtained in the experiments on the Relativistic Heavy Ion Collider (RHIC) [8]. 

Within the present work we provide the detailed study of the linear perturbations of a wide 
class of black holes and find the quasi-normal spectrum for scalar, Dirac, Maxwell fields and the 
gravitational perturbations. We have studied quasi-normal spectra of these fields in four and higher 
dimensions, within the Einstein, Einstein-Aether theories, brane-world scenarios with and without 
the Gauss-Bonnet corrections. We have considered the infiuence of the cosmological constant on 
the quasi-normal frequencies and on the asymptotic behavior of high overtones. We have calculated 
quasi-normal modes of charged scalar and Dirac fields in the Kerr-Newman-de Sitter background. 

Also we have studied quasi-normal spectrum of the massive scalar field for Schwarzschild, 
Tangherlini and scalar hairy anti-de Sitter black holes and the massive vector field in the Schwar- 
zschild background. We provided comprehensive discussions about properties of the quasi-normal 
spectrum and late-time tails of massive fields. We have considered appearing of the infinitely long- 
living oscillations (quasi-resonances) for particular values of the masses of the black hole and of the 
field. 

We prove the instability of the higher dimensional Reissner-Nordstrom-de Sitter black holes for 
sufficiently large values of the black hole charge and of the cosmological constant. By the numerical 
analysis of the quasi-normal spectrum, we support the stability of the Kaluza-Klein black holes with 
squashed horizons, the instability region of the Gauss-Bonnet black holes and the long wavelength 
instability of black strings. We show the appearance of those instabilities in the time domain. Also 
we prove the stability of the massive scalar field in the background of the Kerr black hole. 

We review the numerical methods used for our analysis: the time-domain simulation, the ap- 
proximation by the Poschl- Teller potential, the WKB approach and the Frobenius method. 

We propose a generalisation of the NoUert improvement of the Frobenius method for higher 
dimensional problems. This improvement provides better convergence of the numerical procedure, 
which is crucial for the case when the imaginary part of the quasi-normal frequency is much larger 
than its real part and for the calculation of the quasi-normal modes of massive fields. 

Also we describe a new numerical tool, which we developed for the calculation of the quasi- 
normal frequencies of a black hole, which metric is not known analytically, but can be found as a 
numerical solution of a set of differential equations. This technique supports the natural expectation 
that the dominant quasi-normal frequencies depend mainly on the region nearby the black hole and 
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do not depend on the behavior of the metric at large distance. We have also checked this method 
by the time-domain integration. 

In addition, we discuss an interesting possibility of observation of the acoustic analogue of the 
Schwarzschild black hole in laboratories. 

This work is aimed to clarify a number of questions about the behavior of perturbations of black 
holes. We provide a comprehensive analysis of the black hole perturbations in various theories in 
order to give the complete picture of quasi-normal modes, late-time tails and stability of black holes. 
The presented results were published inPllQlIIIlinillSliailSlIISlIIZlIIHllISlEniE^ 
[251 121]. 
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Chapter 2 

Black hole perturbations 



2.1 Equations of perturbations near a black hole 

The dynamics of the general relativity in D space-time dimensions is described by the Einstein- 
Hilbert action 

S = [ VWl (^(R - 2A) + Ui] d^'x, (2.1) 



where we use the geometrized unit system, so that the speed of light, c, and the gravitational 

constant, 7, are set equal to one. The metric signature is chosen as (H . . .). is the 

Ricci scalar. Cm describes all matter fields appearing in the theory, and A is the cosmological 
constant. 



The variation of the action (2.1) allows to find the Einstein equations 

Rab — ^Rgab + ^dab = ^T^Tab, 

and the field equations in the curved space-time 

i) 



(2.2) 



(2.3) 



The energy-momentum tensor Tab in (2.2) is defined as 

^ 2 d^/^Lu 



- ab 



The general solution of the equations (2.2, 2.3) is some metric and fields in the background this 
metric 

ds'^ = gab{x)dx''dx'', = 4>{i)ix). (2.4) 

We study a stationary black hole, for which gab and do not depend on time in the appropriate 
coordinate system and gauge. 

Let us consider perturbations of the metric and fields as a sum of the unperturbed background 
g'^fj, 0°-) and the actual perturbations Sgat, ^(p^i) 



9ab 



9ab + 8 gab] 
AO 



(i) +0</>(i)- 



(2.5) 
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The perturbations are assumed to be small, i. e. we neglect the contributions of order 0{6gab)'^, 
0{5gab5(t){i)), 0{5(t){i)Y and higher. 

2.3) and taking into account that g^^^ and 0°^^ satisfy the 
e to find the set of hnear equations for the perturbations 5gab 



After substituting (2.5) into (2.2 



equations (2.2, 2.3) as well, we are ab 
and 

In this chapter we consider the stationary spherically symmetric solution of the equation (2.2) 
in 4 dimensions with Tab = 0, which describes a black hole in vacuum. This solution is given by the 
well-known metric 



ds^ = fir)de - -7^- - r\de^ + sin^ 



fir) = 1 A- 

r 3 



(2.6) 



where M is the black hole mass and A is the cosmological constant. 

For A = the metric (2.6) describes a Schwarzschild black hole. If the cosmological constant 
is positive (A > 0), we obtain a Schwarzschild-de Sitter black hole, while for negative values of the 
cosmological constant (A < 0) the black hole is called Schwarzschild-anti-de Sitter. 

Since Tab = 0, the field perturbations in such background are not coupled to the perturbations 
of the metric and, therefore, are equivalent to the test fields in the black hole background. 

Let us describe the variable separation. As an example we consider the test scalar field which 
satisfies the Klein-Gordon equation in the curved space-time 



d 



ab 



- d<^{x) _ _J_d^ _ d_^^^^d^ _ 2/(r) _ /\e,^^ _ ^ 



(2.7) 



where 



A, 



I d . ,d 1 92 
smfc*— + 



is the angular part of the Laplasian. 

In order to separate the angular variables we choose the following ansatz: 



1=0 m=-l 

where Yi^m[6, 0) are the spherical harmonics, which are eigenf unctions of the operator Ag^, 

A,,^il,„(^,0) = -/(/ + i)yi,„(^,0). 



(2.8) 



(2.9) 



The integers / > and \m\ < I are called the multipole number and the azimuthal number respec- 
tively. 



After substitution (2.8) into (2.7) we find the wave-like equation for the function '^i^rnit,r): 



dt 



2 ^^.+Viir)]^Ut,r) = 0, 



(2.10) 



where the effective potential has the form 

Vi{r) = fir 



/ = 0,1,2 



The variable r^, is called the tortoise coordinate. It is defined up to an arbitrary constant as 

dr 



(2.11) 



dr^ 



fir)' 



(2.12) 
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If A > the tortoise coordinate maps all the region outside the black hole to the interval 
[—00, 00). In the anti-de Sitter space (A < 0) the interval becomes (— cx),r^oo), where r^^oo is a finite 



value that depends on the integration constant of (2.12). 



We can also consider other test fields near the black hole. In the most cases, after separation of 



the angular variables, we are able to reduce the equations of motion for such fields (2.3) to the set 



of the wave-like equations of the form (2.10). 



One could find that the test massless Dirac field is described by the effective potentials [29] 



V, 



D± 



fir) 



K 



1,2,3... 



and the electromagnetic field is described by the potential 

/(/ + 1) 



EM 



fir) 



I 



1,2,3 



(2.13) 



(2.14) 



(2.5) 



The equation (2.2) allows us to find the equations of motion for the metric perturbations Sgab 
It is convenient to classify the tensor components 6gab with respect to the transformation 
law under rotations on the sphere around the black hole. They can be of scalar, vector and tensor 
type, denoted by s, v and t respectively 



5g, 



ab 



/ 


s 


s 


V 


V 


\ 




s 


s 


V 


V 






V 


V 


t 


t 




\ 


V 


V 


t 


t 


) 



(2.15) 



In order to simplify the equations for the metric perturbations we can use the invariance under 
infinitesimal coordinate transformations. These transformations act as gauge transformations upon 
the metric perturbations 



X 



X 



+ 5x' 



Sgab Sgab + SXa-b + SXb-a- 



The gauge freedom allows to simplify perturbation equations [31j- After the Regge- Wheeler gauge 
fixing, the perturbations are described by the potentials [30] 



Vsir) 
VJr) 



2/(r) 9M3 + 3c^Mr^ + c\l + c)r^ + 3M\3cr - Ar^ 



(3M + cr)2 



fir) 



lil + 1) 



QM 



lil + 1] 



(2.16) 
(2.17) 



/ = 2,3,4... 



for the scalar and vector types respectively. The perturbations of the tensor type in 4 dimensions 
are eliminated by this gauge fixing. However, the perturbations of the higher dimensional black 
hole metric can be of all the three types [32] . 

According to the Chandrasekhar [33] classification, the perturbations of vector type are called 
axial because they impart the differential rotation to the black hole. Perturbations of scalar type 
are called polar. They are related with infinitesimal deformations of the event horizon. Any pertur- 
bations of a spherically symmetric black hole in 4 dimensions can be divided to its polar and axial 
part in any gauge 
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One should note that the multipole numbers / = 0, 1 for gravitational perturbations and / = 
for the electromagnetic field perturbations are not dynamical. For instance, / = gravitational 
perturbations are spherically symmetric and, therefore, obey the Birkhoff theorem. They correspond 
to infinitesimal change of the black hole mass. Similarly, / = perturbations of the electromagnetic 
field and / = 1 gravitational perturbations correspond, respectively, to the infinitesimal shift of the 
charge and the position of the black hole pTj . 



2.2 Time evolution of perturbations 



In order to integrate the equation ( 2.10[ ) numerically we use the technique developed by Gundlach, 
Price and PuUin |35]. We rewrite the wave-like equation (2.10) in terms of the so called light-cone 
coordinates du = dt — dr^, and dv = dt + dr^. 



dudv 



+ V{u,v) I ^(M,f) = 0. 



Let us consider the operator of time evolution in these coordinates 



exp ( h 



d_ 



exp h— + h — 
ou ov 



exp(/.|-)+exp(;.|- 



dudv 



Acting by this operator on \E', and taking into account (2.18), we find 



0{h^). 



(2.18) 



(2.19) 



where we introduced letters to mark the points as S 
N = {u + h,v + h). 



u,v) 



W 



{u + h,v), E = {u,v + h), 



The equation (2.19) allows us to calculate the values of \1/ inside the rhombus, which is built on 
the two null-surfaces u = uq and v = vq (see fig. 2.1), starting from the initial data specified on 
them. 

Let us study, as a qualitative example, the time-domain profile of the vector type gravitational 
perturbations of the Schwarzschild black hole \l/(t,r = R) (fig. 2.2). We can divide the evolution 
of the perturbations into three stages. The first stage depends on the initial conditions and on the 
point R. At late time {t ~ R) we see exponential damping of the amplitude of the perturbations, 
which is followed by the so-called tails at asymptotically late time (t ^ M). 

One can observe, that, after the initial outburst at the first stage, the behavior of the amplitude 
of the perturbations does not depend on the initial conditions. Being independent on the source of 
the perturbations, it depends only on the parameters of the field and the black hole. Therefore, the 
late-time damping law appears to be an important characteristic of the black hole. 

The exponential damping of the perturbations is called quasi-normal ringing. It can be split 
to the superposition of exponentially damping oscillations, that can be represented as a set of 
complex frequencies, which are called quasi-normal modes. The real part of a quasi-normal frequency 
describes the actual frequency of the oscillation, while the imaginary part is its damping rate. In 
the next section we will study the basic properties of the quasi-normal spectrum of black holes. 

In the case under consideration, the late-time tails decay according to the inverse power law. 
The qualitatively similar behavior is observed for the decay of all massless fields in asymptotically 
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Figure 2.1: The integration grid. Each cell of the grid represents an integration step. The thick 
points illustrate the choice of (5, W , E, N) for the particular step of the integration. The initial 
data are specified on the left and bottom sides of the rhombus. 



fiat backgrounds. Yet, the late-time tails appear to be very sensitive to the asymptotical behavior 
of the potential: in the asymptotically anti-de Sitter space-times the quasi-normal ringing governs 
the decay of perturbations at all times [7] . The black hole perturbations in the de Sitter space-time 
have exponential tails, if the cosmological constant is large, and both power and exponential tails, 
if the cosmological constant is small [361 137] . 
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Figure 2.2: The three stages of the evolution of the Schwarzschild black hole gravitational pertur- 
bations {1 = 2 vector type, A = 0, R = 22M). Time is measured in units of the horizon radius. 



2.3 Properties of the quasi-normal spectrum 
2.3.1 Boundary conditions 

In order to study quasi-normal spectrum of a black hole it is convenient to make Fourier transform 
for the function 



-iuJnt)Qu 



(2.20) 



n=0 

Let us note that generally the functions {exp(— 2ti;„)}, where {ujn} is the quasi-normal spectrum, 



do not form basis in the vector space of the solutions of the equation (2.10). Therefore, the signal 



cannot be expanded in terms of these functions for all time. Indeed, as we see on fig. |2.2[ the 
exponential decay is replaced by slower power-law decay at asymptotically late time. Yet, the 



expansion (2.20) is appropriate to describe the quasi-normal ringing epoch. 



Henceforward, we will omit the index n of cj, implying that to is any frequency from the quasi- 
normal spectrum. 

The function satisfies the linear equation 



+ u:'-V{r)\Q^{r) = Q. 



(2.21) 



Thus, the searching of the quasi-normal modes is reduced to the problem of finding eigenfre- 



quencies of the equation (2.21 ). Since we are interested in what happens in the i?-region outside the 



20 



black hole, we have to impose the boundary conditions at the event horizon and the cosmological 
horizon for A > (or the spatial infinity for A < 0). In order do this we must study the structure 
of the singularities of the equation (2.21) at those points. 



Let us study the general properties of the singular points of the equation (2.21). First we note 
that at any horizon r ^ h 

V oc /(r) =o(-]. (2.22) 



To prove (|2.22|) we consider two types of the horizon: 

-1 (l + o(r-/i)). 



1. /(r -,h) = {r- h)f'{h) + o(r - h) 
then ^ 

= -irTTT In 



h 



Thus, we find that the dominant contribution is 

/(r ^ h) hf'Qi) exp {f'{h)n) ~ hf'{h) exp ^In 

This term decays exponentially with respect to because 

- 1 



lim In 

i — >h 



— OO. 



2. /(r ^h) = A{r - hY + o(r - /i)", a > 1 then 

1 



and 



/(r ^ h) 



A{a - l)(r - h) 
r — h 



-{l + o{r-h)) 



(a - l)n 



+ o(r-/i)° = 



If A = 0, the condition (2.22) at spatial infinity is also satisfied for the potentials (2.11, 2.13 



2.14| |2.16[ |2.17D . Indeed, r ~ and 

V{r) = O 



O 



Since the potential satisfies (2.22), the behavior of the eigenfrequency at the boundaries can be 

oc exp(±iu;r^). (2.23) 

The appropriate boundary conditions for the problem under consideration are purely ingoing wave 
at the event horizon and outgoing wave at the cosmological horizon (spatial infinity) 



-OO : 
-OO : 



oc exp(-ico'r^), 
oc exp{+iujr^). 



(2.24) 



This choice of the boundary conditions imply that, after the black hole is perturbed, there is no 
signal, which comes from the black hole or from any remote source [33]. It means that the source 
of the perturbations acts only before the first stage of the evolution. After this, we consider only 
the black hole response upon the perturbations. 
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For the 4-dimensional Schwarzschild-anti-de Sitter black hole (A < 0), one can find that the 
linear independent solutions at spatial infinity (r oo) are 

Qi oc r, Q20C r~^. 

Since we do not suppose to have an infinite amplitude of the perturbation at the spatial infinity, 
we must choose Q2 as the appropriate boundary condition there, i. e. we impose the Dirichlet 
condition at the AdS boundary [7J 

QUr = 00) = 0. (2.25) 

The Dirichlet boundary conditions are usually imposed at the spatial infinity in asymptotically 
anti-de Sitter backgrounds (see e. g. [38j). 



2.3.2 Quasi-normal spectrum and black hole stability 



Black holes do not exist in nature, that is why, the important property of any black hole solution 
is its stability against perturbations. Unfortunately, it is possible to prove stability analytically 
only of some relatively simple solutions of the Einstein equations. That is why the numerical test 
of stability is important for the black hole study. The instability implies the existing of growing 
modes in the quasi-normal spectrum. It is clear, that the linear approximation is enough for this 
test. Indeed, if the black hole is stable, any perturbations will decay until the linear approximation 
is valid. If the quasi-normal spectrum has a growing mode, the amplitude of perturbations will 
grow until we are compelled to consider the non-linear back reaction of the perturbations upon the 
metric. In this case, we state, at least, that in order to obtain the stable solution we must take into 
account this non-linear correction, which could be also non-stationary. 

It is important to note that the black hole is unstable if there is only one growing mode in its 
spectrum. Therefore, in order to prove the black hole stability we must show that the quasi-normal 
spectrum does not contain any growing mode, i. e. frequency with positive imaginary part. It 
makes the numerical proof of the black hole stability extremely complicated. 

Yet, in some cases we are able to prove the stability analytically. Let us multiply the equation 



(2.21 ) by the complex conjugated and integrate the first term by parts 




co'\QUn)\'-V\QUr.)\'- 



dQu 



dr^ 



dn = 0. 



Taking into account the boundary conditions (|2.24|) we find 
iuA + u^B = 




Vir)\Q^in)\^ 



dQuj{n 



dr^ 



dr^ 



(2.26) 



00 

where A = |g^(r, = oo)\^ + \Q^in = -00) ^ > 0, B = j \Q^{n)\^dr, > 0. 



The imaginary part of (2.26) reads 

Re{uj)A + 2Re{uj)lm{uj)B = 0. 



(2.27) 



We can see that lm{uj) < except for the case Re(c(j) = 0. Thus, we conclude that the growing 



modes do not oscillate (see sec. 5.2). 
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If the potential is positive in the i?-region outside the event horizon and Re(ct;) = 0, the real 
part of (2.26) implies that the imaginary part of the frequency remains negative 

-lm{uj)A - lm{ujfB > Im(tu) < 0. 



It is clear, that the righthand side of (2.26) can be positive even though V{r) is not positive 
everywhere. By introducing the new derivative D = h S{r^) we can rewrite the integral as 



Vir)\QUn)\' + 



dQu 



dr^ 



V\QUr.)\^ + \DQ^ndr 



S{n)\QUn)\' 



where V 



dr^ 



Thus, we conclude that if we find the function S'(r^) such as S{r^, = oo) < 0, 



S{r^, = — oo) > and > 0, the righthand side of (2.26) stays positive and, therefore, Im(ci;) < 0. 

This technique is called S- deformation and allows us to prove stability for some cases, when the 
potential is not positive definite [39] . It is important to note that if we find an appropriate function 
S, we prove the black hole stability. Otherwise we do not know if the black hole is stable or not. 

The same approach can be used in a similar way for anti-de Sitter black holes. We must just put 
r^oo as the upper bound in the integrals and use (2.25). Then we obtain A = \Quj{r^ = — C)o)p > 
and the requirement for 5* at the spatial infinity reads S{r^ = r^ao) < oo. 

Because we must guess the appropriate function S, the technique of S-deformation allows to 
proof stability only for relatively simple potentials. For more complicated cases we are able to prove 
the black hole stability only numerically. 



2.3.3 Isospectrality 



Let us consider two equations (2.21) with the effective potentials, taken as and V 

dW{r, 



with 



± 



dr^. 



(2.28) 

where H^(r^) is some finite function and /3 is a constant. Then, if is an eigenfunction of (2.21) 
for the potential , the eigenfunction for the potential V~ is given (up to an arbitrary factor) by 

d 

dr^, 



q: 



w 



(2.29) 



corresponding to the same eigenvalue uj. Thus we conclude that the quasi-normal spectrum is the 
same for the potentials and V~ . 

Let us consider the examples of isospectrality for 4-dimensional Schwarzschild((-anti)-de Sitter) 
black holes. One can check [33j that the potentials (2.16) and (2.17) can be obtained by taking 



2M 3 + 2c 3c2 + 2c3-9AM2 
W = —^ ^ + 



3M 



& + c- 



3AM2 



/5 



If 



9M2 



3c(3M + cr) 

The potentials for the Dirac field ( |2.13[ ) are, obviously, isospectral too. 

It is clear that if one of the isospectral potentials is positive outside the event horizon, as one 
can notice for Vs{r) (c > 2), both do not lead to growing modes, implying stability against the 
perturbations. The other possible case of stable potentials is /3 > 0, as it happens for (2.13). In 
order to prove the stability we just use the S'-deformation technique with S = ±W. In fact, all the 
potentials (2.11, 2.13, 2.14, 2.16, 2.17) are positive for — oo < < oo, thereby, the stability of the 
Schwarzschild black hole is evident. 
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Chapter 3 



Numerical methods of calculation of the 
quasi-normal modes 



3.1 Fitting time-domain data 

The most direct approach to finding quasi-normal modes is the numerical integration of the equation 
(2.10) as described in the section 2.2 The result of the time-domain integration is a time profile 
data {^'(t = 0), \E'(t = h),'^(t = 2h) . . .}, which can be used to calculate the quasi- normal modes. 

Let us describe the simplest Prony method of fitting the profile data by superposition of damping 
exponents (see e. g. |10] and references therein) 

p 

^ ^ae-*^»*. (3.1) 

i=l 

We suppose that the quasi-normal ringing epoch starts at t = and ends at t = Nh, where 
integer > 2p — 1. Then the formula (3.1) is valid for each value from the profile data 

p p 

The Prony method allows to find Zi in terms of known x„ and, since h is also known, to calculate 
the quasi-normal frequencies cuj. In order to do this, we define a polynomial function A[z) as 

p p 

A{z) = l[{z-z,) = J2(^n.zP-^, ao = l. (3.3) 

j=l m=0 

Let us consider the sum 

p p p p p p 

J2 = E «™ E c^^r"" = E c^^r E «-^r"^ = E ^.^r'^(^.) = q- 

m=0 m=0 j=l j=l m=0 j=l 

Taking into account that ao = 1, we find 

p 



m=l 
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Substituting n = p..N into (3.4) we obtain N — p+l > p linear equations for p unknown coefficients 
Let us rewrite these equations in the matrix form 



Xp—i Xp—2 
Xp Xp—i 



Xo 
X\ 



\ /«1 \ 

Oi2 



y XAr_i XAr_2 . . . ^AT-p j \ j 



Xp 



Xp^\ 



\ Xn J 



Such matrix equation 



Xa = —X 



can be solved in the least-squares sense 



(3.5) 



where denotes the Hermitian transposition of the matrix X. 

After the coefficients of the polynomial function A{z) are found, we can calculate numerically 
the roots Zj of the polynom and the quasi-normal frequencies 



h 



Because the quasi-normal stage is not a precisely defined time interval, in practice, it is difficult 
to determine when the quasi-normal ringing begins. In fact, when we observe explicitly damped 
oscillations, we usually see only the fundamental mode, while higher overtones, which damp quickly, 
are already exponentially suppressed. Being a small corrections to the signal, such higher damped 
oscillations are indistinguishable from numerical errors within the described approach. Thus, the 
higher overtones are difficult to find. Usually, the Prony method allows to calculate at most six roots 
of the polynom A{z), including the complex conjugated ones, which correspond to the symmetry 
{—uj*)*. This symmetry exists just because the wave-like equation (2.10) is real. Therefore, in 



UJ 



fact, we are able to calculate only two or, sometimes, three dominant frequencies. 

In order to determine the beginning of the quasi-normal ringing epoch more precisely, we can use 
the following technique [H] . Let us find the dominant quasi-normal mode uJi and the corresponding 
coefficient Ci at some late time interval. Then we can subtract this oscillation from the numerical 
data, and obtain the proffie data without the contribution of the dominant mode. After this, one 
can see ringing for the ffist overtone. If the lifetime of the quasi-normal ringing is long enough, we 
are able to find, step by step, higher overtones, making sure, that the numerical error of the initial 
data is less than the signal after removing the contributions of the lower-damping modes. 



3.2 Approximation by the Poschl- Teller potential 

The easiest method of calculation of the quasi-normal modes in frequency domain is approximation 
of the effective potential by the Poschl- Teller potential. This method was suggested by Bahram 
Mashhoon \M 



Suppose that the potential in the equation (2.21) is invariant under the following transformation 



V{r\a) = V{-ir*,a'), 
where a is some parameter and a' depends on a. 
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Let us consider the solution of the equation with the invese potential 



^ + i-Q' + V)Qn = 0, (3.6) 

with the boundary conditions that are characteristic of bound states if Re(fi) > 

Qnoce^^'-*, n^±oo. (3.7) 



It is easy to see, that this solution is related with the solution of (2.21) with the quasi-normal 



boundary conditions (2.24) in a simple way 

QUr*,a) = Qn{-'ir*,a'), ^{a) = Q{a'). (3.8) 

Thus, the quasi-normal mode problem is reduced now to the bound states problem for an inverse 
potential V —V, which is smooth potential gap, approaching zero at the infinite boundaries. This 
gap can be approximated by the Poschl- Teller potential 

VpT = ^, (3.9) 

cosh a(r=K — rx) 

where Vq is the height of the effective potential and — 2Voa^ is the curvature of the potential at its 
maximum. The bound states of the Poschl- Teller potential are known in analytical form [44j 

n = «|-[n + i] + Q + ^] ^ 1, n = 0,l,2,... (3.10) 



The potential (3.9) is obviously invariant under the transformation 

VpT{r^,a) = VpT{-ir^,ia). 
Therefore, the quasi-normal modes for the inverted Poschl- Teller potential are 

2\ 1/2 / 1 

a \ . / 1 



uj{a) = n{ia) = ± - — j -ia\^n+-j, n = 0,1,2,... (3.11) 

Technically one has to fit a given effective potential to an inverted Poschl- Teller potential. In 
the chapter [7] one shall see that in many cases the behavior of the effective potential only near 
the black hole is essential for determining the dominating quasi-normal modes. So that the fit of 
the height of the effective potential Vq and of its curvature — 2Vo«^ is indeed sufficient to estimate 
quasi-normal frequencies. 

This method gives quite accurate results for the regime of high multipole numbers I. In par- 
ticular, for gravitational perturbations of the D = 4 Schwarzschild black holes, the fundamental 
quasi-normal modes obtained by the Mashhoon method gives relative error of not more than 2% 
for the lowest multi-pole I = 2, and of about fractions of one percent for higher multipoles. 

There are cases when the effective potential of a black hole is exactly the Poschl- Teller potential. 
These are Schwarzschild- de Sitter |15] and Reissner-Nordstrom-de Sitter black holes with extremal 
value of the A-term in D > 4 space-time dimensions 
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Figure 3.1: The three regions separated by the two turning points Q{x) = 0. 

3.3 WKB method 



In order to evaluate quasi-normal modes for more complicated effective potentials it is convenient 
to use the WKB (Wentzel, Kramers, Brillouin) method, which provides good accuracy. The WKB 
technique was applied to finding of the quasi-normal modes of black holes for the first time by 
Schutz and Will [17]. 



In order to simplify our notations, let us re- write the wave-like equation (2.21) in the following 
form ^ 

^ + Q{x)^ix) = 0, (3.12) 

i.e. we identify x = r^:, Q ^ u"^ — V, and = 

Let us introduce the WKB parameter e in order to track orders of the WKB expansion. The 
asymptotic WKB approximation at both infinities has the following general form 



oc exp j 



(3.13) 



vn=0 



Substituting the expansion (3.13) into the wave equation (3.12), and equating the same powers of 
e, we find 



So{x) = ±z j Q{7^y/^dv, 

XQ 



Si(x) 



InQ(x). 



(3.14) 



(3.15) 



The two choices of the sign in (3.14) correspond to either incoming or outgoing waves at either of 
the infinities x = ±oo. 

Thus, at X — s> +00 (region I), Q{x) — > cj^ in the dominant order, so that Sq +iu!x for a wave 
outgoing at the infinity and 5*0 —iujx for a wave in-coming from infinity. In a similar fashion, at 
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the event horizon x —>■ —oo (region III), 5*0 +iujx is for wave in-coming from the event horizon, 
while 5*0 —iujx is for a wave out-going to the event horizon. We shall designate these four 
solutions as \E':'^, \E'i, '^^^^ and \E'{^^ respectively for plus and minus signs in 5*0 in I and III regions 
(see fig. 3.1 ). Thus 

^/ _ ^+iu>x^ ^ _^ ^7 _ ^-iux^ x^+oo (3.16) 

^in^^+iux^ a;^-oo, ^^^^ ~ e"'"", x -> -oo (3.17) 
The general solutions in the regions I and III are 

m = Zlj>i + Zi^,¥^, (3.18) 
* = C'^'/' + C'/^- (3.19) 

The amplitudes at a; — >■ -|-oo are connected with the amplitudes at x — »• — oo through the linear 
matrix relation 

(3.20) 



^in J \ ^21 S22 J V ^in 



Now we need to match both WKB solutions of the form (3.13) in the regions I and III with a 
solution in region II, through the two turning points Q{x) = 0. 

If the turning points are closely spaced, i.e. if —Q{x)max ^ ^(ioo), then the solution in the 
region II can be well approximated by the Taylor series 

Qi^) = Qo + ^Qoix - a;o)' + O ((x - x^f) , (3.21) 

where xq is the point of maximum of the function Q{x\ Qq = Q{xo), and Qq is the second derivative 
with respect to x taken at the point x = xq. Region II corresponds to 



\x-Xo\<,l^^ey^ (3.22) 

The latter relation gives also the region of validity of the WKB approximation: e must be small. 
Let us introduce new functions 

t = {2Q';^'/'e'-/\x-xo), (3.23) 
^+1 = -iQo/i'iQlf". (3.24) 



Then the wave equation (3.12) takes the form 

(fm f 11 



+ z/ + ---r^ = 0. (3.25) 



The general solution of this equation can be expressed in terms of parabolic cylinder functions 

m = AD^{t) + BD^^^iiit). (3.26) 
Large \t\ asymptotics of this solution are 

^ ^ 5e-3''-('^+i)/^(4A:)-('^+i)/^(x -xo)-('^+^)e^'='''(^-^°)'/2^ 
{A + B(27r)i/2e— /2^^^^ ^ l))e''^'^/^(4A;)'^/^(a; - a;o)"e-*'='^'("-"»)'/^ x > X2, (3.27) 
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{B - iA{27T 
1 



,1/2, 



-W2/r(_^))gi-(-+i)/4(4^^ 



X — Xo 



-ik^/^(x-xo)^/2 



+ 



X - Xo) 



-xof/2 



X < xi, (3.28) 



where A; 



Equating the corresponding coefficients in (3.27), 
the elements of the S matrix, 



(3.28) and ehminating A and B, we obtain 



yiii 

^out 
ylll 



R-\2Txf/^/T{-v) 



2i?V^'^(27r)i/Vr(z/ + l) 



7^ 

^out 



where 



R 



(;.+l/2)/2^-(i/+l/2)/2 



(3.29) 



(3.30) 



When expanding in higher WKB orders, the 5* matrix has the same general form with just other 
i?, still depending only on v. Let us note that for a black hole there are no waves "reflected by the 
horizon", so that Z^^^ = 0, and due to quasi-normal mode boundary conditions, there are no waves 
coming from infinity, i.e. Z- 



0. Both these conditions are satisfied by (3.29), only if 
Ti-u) 



oo, 



and, consequently, u must be integer. Then, from the relation (3.24) we find 

1 



n + 



-zQo/{2Q'^'/\ 



n 



0,1,2, 



(3.31) 



(3.32) 



The latter relation gives us the complex quasi-normal modes labeled by an overtone number n with 
the accuracy of the ffist WKB order [17]. Later this approach was extended to the third WKB 
order beyond the eikonal approximation by Iyer and Will [38] and to the sixth order by Konoplya 
M l50] . In order to make the higher order WKB extension it is sufficient to take higher orders in 



e WKB series (3.13) and to take appropriate number of consequent terms in the Taylor expansion 
(3.21 ). Since the S-matrix (3.20) depends only on z/, its elements Sij can be found simply by solving 
the interior problem in region II at higher orders in e [IB], and without explicit matching of the 
interior solution with WKB solutions in regions I and III in each order. 

Going over from Q to the effective potential V , the sixth order WKB formula has the form 

'^""IZ^^ - A2 - A3 - A4 - As - Ae = ^ + ^, n = 0, 1, 2, .. (3.33) 
a/-21/o / 

where the correction terms Aj depend on the value of the effective potential and its derivatives (up 
to the 2i-th order) in the maximum. 

It was shown in [51] that WKB formula, extended to the sixth order, gives the relative error 
which is about two orders less than that of the 3rd WKB order. Yet one should remember that 
strictly speaking the WKB series converge only asymptotically and the consequent decreasing of the 
relative error in each WKB order is not guaranteed. Therefore it is reasonable to develop a modified 
WKB technique in the so-called optimal order [32]. The latter gives better results for moderately 
higher overtones n and especially when n > I. Yet, in many cases when n < I the 6th order WKB 
formula gives better results than the optimal order treatment. 

In some cases, the WKB approach needs modifications: for instance when considering a massive 
scalar field in a black hole background, the effective potential has a local minimum far from a black 
hole. This local minimum induces two changes in the WKB procedure. First, there are three turning 
points which separate all space into four regions, so that three matchings are required. Second, an 
influent subdominant term in the asymptotic WKB expansion at spatial infinity (3.19) appears (see 
sec. IsXtI). 
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3.4 Frobenius method 



3.4.1 Frobenius series 



The most accurate method of searching of eigenvalues of the equation (2.21) is Frobenius method, 



which allows to find quasi-normal modes with arbitrary precision and does not require special form 
of the effective potential. This method allows to calculate also higher overtones of black holes. It 
was done for the first time by E. W. Leaver for Schwarzschild and Kerr black holes 



(3.34) 



Let us consider the second order differential equation of more general form than (2.10) 



+ p{r)— + q{r) ) R{r) = 



where the functions p{r) and q{r) depend on the eigenfrequency u. 

Let us start from the analysis of character of singular points of this equation. There are two 
points, which are always singular: the event horizon r = r_|_ and the cosmological horizon (or the 
spatial infinity) r = r^o- Usually, there are also other singular points, that depend on p{r) and 
q{r). By definition, the quasi-normal modes are eigenvalues of uj with the boundary conditions that 
correspond to the outgoing wave at the cosmological horizon (spatial infinity) and the ingoing wave 
at the black hole event horizon. So, we are able to define the function R{r) as a multiplication of 
some factor and the Frobenius series. The factor is divergent at these singular points. It is chosen in 
order to the series be convergent in the region r_|_ < r < Too- If p{r) and q{r) are rational functions 
of r, we can construct such series in terms of the rational functions, 



R{r) 



in 



r — r_i 



r — Ta 



r 



-if} oo 

fc=0 

oo 

fc=0 



r 



r - ro Too 

k 

r — 



r 



OO. 



(3.35) 



and 



The values f2, a and /? are defined in order to satisfy (3.34) in the singular points r = 
r = roo- The quasi-normal boundary conditions fix Re{Q) and Re{(3), which must be chosen of the 
same sign as Re{u). 

Let us consider the series 

oo 



(3.36) 



k=0 



If all the singular points of the equation ( 3.34^ satisfy \z\ > 1, the series (3.36) are convergent at 
z = 1 [r = Too), if and only if the value of uj is the eigenfrequency of the equation (3.34). If there 



is at least one singular point inside the unit circle, one has to continue the series (3.36) through 



some midpoints (see sec. 3.4.4) in order to test the convergence at the cosmological horizon or at 
the spatial infinity. 

Note, that the definition of z contains an arbitrary parameter ro < r+. In most cases, it can be 
chosen in order to move all the singularities outside the unit circle. 



3.4.2 Method of a continued fraction 



Substituting (3.35) into (3.34), one can obtain an A^-term recurrence relation for the coefficients bi 

min{N—l,i) 

-'^^\uj)bi.j = 0, i>0, (3.37) 



E 

j=0 



-'3,1 
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where the coefficients CjJ\uj) (0 < j < min{N — l,i)) depend on uj. 
We now decrease the number of terms in the recurrence relation 



min{k,i) 
3=0 



(3.38) 



by one, i. e. we find cf^{uj), which satisfy the equation 



min{k—l,i) 

i=o 

For i > k, we can rewrite the above expression as 



0. 



(3.39) 



(3.40) 



Subtracting (|3.40|) from (|3.38|) we find the relation (|3.39|) explicitly. Thus we obtain, 

0, or i < k. 



J 



C,,i [UJ)Cj_,i_^[UJ) 



This procedure is called Gaussian eliminations, and allows us to determine the coefficients in the 
three-term recurrence relation numerically for a given lj up to any finite i 

cS b^ + 45 h-1 + 45 = 0, ^ > 1 (3.41a) 
45 bi + 45 bo = 0. (3.41b) 



The complexity of the procedure is linear with respect to i and A^. 



If the the Frobenius series are convergent, we are able to find 61/60 from (3.41b) and substitute 

(3.42) 



it into (3.41a) 



bo 



J3) 
"0,1 



^(3) „(3)^(3) (3) (3) 
•-2,2 '-0,2'-2,3 '-0,3'-2,4 



.(3) 



.(3) 



.(3) 



'^1,2 '^1,3 '^lA 



Finally we find 







'"1,1 



„{3) J3) „(3) (3) 
'-0,l'-2,2 '-0,2'-2,3 



.(3) 
-1,2' 



.(3) 
"1,3' 



what can be inverted n times 

-.(3) J3) 



'"l,n+l 



.(3) ^(3) 



-.(3) 

-'l,n-l' 



.(3) 

-l,n-2' 



„(3)„(3) 
'-2,2'-0,l 

'^1,1 



„(3) „(3) (3) (3) 
'-0,n+l'-2,n+2 '-0,n+2'-2,n+3 



.(3) 

--1,11+2' 



.(3) 

-l,n+3' 



(3.43) 



(3.44) 



The equation (3.44) with the infinite continued fraction on the right-hand side can be solved 
numerically by minimising the absolute value of the difference between the left- and right-hand 
sides. The equation has an infinite number of roots (corresponding to the quasi-normal spectrum), 
but for each n, the most stable root is different. In general, we have to use the n times inverted 
equation to find the n-th quasi-normal mode. The requirement that the continued fraction be itself 
convergent allows us to restrict its depth by some large value, always ensuring that an increase in 
this value does not change the final results within the desired precision. 
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3.4.3 Nollert improvement 

It turns out, that the convergence of the infinite continued fraction becomes worse, if the imaginary 
part of uj increases with respect to the real part. It means that in order to calculate the higher over- 
tones correctly, we must increase the depth of the continued fraction, what dramatically increases 
the time of calculation. The convergence is bad also if rg in (3.35) is not a singular point. Such 
fixing of ro is necessary to move all the singular points outside the unit circle for higher- dimensional 
Schwarzschild black holes. 

The problem of slow convergence was circumvented in |54j for the three-term recurrence relation 
and generalised for higher N in [18]. Let us consider 



hn-l 

that for large n can be expanded as 

RJuj) -- 



Rn 



M J3) (3) 



.(3) 



.(3) 



n 



n 



(3.45) 



(3.46) 



In order to find the coefficients Cj of (3.46), we divide the equation (3.37) by hi-N+i and use 
the definition = — &„/6„_i. We find the equation with respect to Rn 



j=0 



N-2-j 

n « 

k=0 



i—k 



0. 



(3.47) 



For large n, c^^^{uo) oc n^. Thus, substituting the expansion (3.46) into (3.47), we find 



lim 



N-l-j. 



0. 



(3.48) 



j=0 



In general, the equation (3.48) has — 1 roots (in fact there are multiple roots). One of the 
roots (also multiple) is always Cq = —1, implying the unit radius of convergence of the series (3.35). 
Other roots appear due to the existing of additional singular points of the equation (3.34). Thus 
we choose Co = —1. 

After fixing Co = — 1 one can find an equation with respect to Cf. In order to fix the sign of Ci 
we can use the convergence of the series (3.35) at ^ = 1. Therefore, 

lim b„ = 0, 



i. e.$N: \fn > N > \h 



n-l 



Since for large n we have 



-Rn ~ —Co — 



bn-1 

we find out that the real part of Ci cannot be negative. 



n 



After the sign of Ci is fixed, the other coefficients in (3.46) can be found step by step from (3.47) 
without encountering indeterminations. 

Since we can calculate the coefficients Cj, the expansion (3.46) could be used as an initial 
approximation for the "remaining" infinite continued fraction. In order to ensure the convergence 
of (3.46) for a given value of u, one has to start from the found approximation deeply enough inside 
the continued fraction (3.44). The expansion gives a good approximation for Rn. Therefore, the 
required depth is less than it would be, if we started from some arbitrary value. 
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3.4.4 Continuation of the Frobenius series through midpoints 



Let us consider the case when we are unable to fix the parameter tq in (3.35) in such a way that all 
the singularities, except r = r+ and r = Too, move outside the unit circle. In this case there is at 
least one singularity, for which \z\ < 1. This singularity implies smaller radius of convergence for 



the series (3.36). In order to test that the function u{z) is convergent at 2 = 1, we must continue 
the series analytically, by constructing iteratively the expansions of u{z) at some midpoints |55] . 
Namely, we equate the series expansion at two points 



uiz 



(3.49) 



n=0 



n=0 



where ^ = is a midpoint inside the radius of convergence of (3.36) 



The coefficients 6„ also satisfy the A^-term recurrence relation, which could be reduced to the 
three-term one 



g(3)- 



hi + cf] + 47 6,„2 = 0, ^ > 1. 



;;(3) 



(3.50) 



In order to find 61/62, we must use the condition at the event horizon by taking into account (3.49) 

00 00 



60 



(3.51) 



n=0 



n=l 



From (3.41b) and (3.41a) we find the coefficients 6„ and substitute them into (3.50). If 2; = 1 is the 



closest singular point to 2; = 2:0, we obtain the equation with respect to uj as 



61 
60 



;;(3) ~(3)~(3) ~;(3)~;(3) 
•-2,2 '-0,2'-2,3 '-0,3'-2,4 



^(3)_ ~(3)_ 
^1,2 '-1,3 



g(3)_ 



(3.52) 



Otherwise one has to repeat the procedure, by constructing the series (3.49) for the next midpoints 
Zi, 2:2, 2:3, . . ., until the cosmological horizon (or spatial infinity) appears to be inside the radius of 
convergence. 



One should note, if the convergence of the continued fraction on the right-hand side of (3.52) is 



slow, one can use the Nollert improvement. Since the radius of convergence of the Frobenius series 



is now less than one (i? < 1), we must choose Cq = —R ^ in (3.46). 



3.4.5 GeneraUsation of the Frobenius series 



The series expansion (3.36) is not necessarily to be done in terms of power of a rational function 



of r. For some purposes, the more convenient choice is expansion in terms of an other full set of 
functions in the appropriate Hilbert space. Here we consider, as an example, a Kerr-Newman-de 
Sitter black hole, that is described by the line element 



where 



dr^ de^ 
a7 ^ a7 



Ae sin^ 6 
(l + a)V 



[adt - (r^ + a^)dLp]'^ + 



;i + a)V 



A, 
A, 



1 - 



a 



1 + a cos^ 9, 



- 2Mr + 
Aa2 



(dt- asin^edipf, (3.53) 



(3.54) 



a 
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M is the black hole mass, Q is the charge, a is the rotation parameter, and A is the cosmological 
constant. 

After separation of the variables, the angular part of the massless (charged) field equation of 
motion can be reduced to 



- 2ax'^+ 



— 1 + ax^){l - x^)— + A - s i - a + ^ ^ 

\dx ax a 

(^s{am - (1 + a)i)x - ^^^^^ - 2m(l + a)i + s^^ - 



l + a 
1 + ax'^ 



(3.55) 



(1 + afra^ 



(1 + a)(s2 + 2smx)\^ 



'(l + ax2)(l-x2) l-x2 l*^^^^ °' 



where A is the separation constant ^ = auj, x = cos{9), s is the field spin and m is the projection of 
the angular momentum of the field onto the axis of the black hole rotation, < s < 2 and m are 
(half) integers. 

The appropriate series for the function S are [56] 



oo 

S{Z) = Z^^iz - iy%Z - Z,)^^{Z - Z^)Y,bnUn{z) 

n=0 



(3.56) 



where 



z 



a — i X + 1 
2 x-v/a — 



2c 



z(l + iy/a) 
4:Ja 



2^/a 



\m — s\ 



A, 



\m + s\ 



1 fl + a I— 

2 \ \ a 



2 ' " 2 ' 

and the expansion is done in terms of the Jacobi polynomials 

«„(.) = F(-n, n + 7; .) = (_i)n r(2n + a;)n! ^(,_^,^_,) _ 

r(n + 7) 

where w = 2y4i + 2^2 + 1 and 7 = 2Ax + 1. 

The coefficients 6„ in ( |3.56[ ) satisfy the three-term recurrence relation (3.41a, 3.41b) with 



,(3) 



n(n + Ai + T g)(n + 2A2) 

\fa 2{2n + 2Ai + 2A2 + l){n + Ai + A2] 

Jn 



(3.57) 



+ 



2{n + Ai + A2){n + A1 + A2-I) 
>-l) (71 + 2^1 + 2^2) 



1 

"4 



A - 2A1A2 - Ai- A2 + 2{m + ST (2Ai + 1))^ 



w? — 



— s 



i {n-l + Ai + A2± s){n - 1 + 2Ai){n - 1 + 2Ai + 2Ar 
2{2n + 2Ai + 2^2 - 3){n - 1 + Ai + A2) 



a 



(3.58) 



(3.59) 



where 



Jn = {n-l){n+2Ai+2A2){Ai~A2) + {Ai+A2±s+l) {{n - l){n + 2Ai + 2A2) + {2Ai + l){Ai + A, 



:))■ 
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Since the series (3.56) must be convergent at z = 1, we can solve numerically the equation (3.44), 
i find, thereby, the separation constant as a function of frequency: , 

In the Reissner-Nordstrom-de Sitter limit (a 0) the equation (3.44) is reduced to ^ = 0. In 



this case the value of A does not depend on cu 

X = (l - s + 1)(/ + s), I = n + Ai + A2> max{\m\, \s\). 

3.4.6 Probenius series for the radial part of the charged field equation 
in the Kerr-Newman-de Sitter background 



The radial part of the massless (charged) field equation of motion is of the form (3.34) and reads 



2a 



4zs(l + a)ujr - —{s + l)(2s + l)r^ + 2s(l -a)- 2isqQ - A }R{r) = 0, (3.60) 



where K = [c<j(r^ + a^) — am](l + a) — qQr, q is the field charge. 
The appropriate Frobenius series are found to be ^ 

R{r) 



r — r. 



r — r^-Vr. 



r — 



(3.61) 



Note, that, in order to obtain the recurrence relation for both types of the boundaries (asymptot- 

dBir) K 

ically fiat and de Sitter), we introduce the exponent e*'^*-'"-' such that — — = — . This exponent 

describes an outgoing wave at the horizons and spatial infinity. Thus we have to compensate the 
outgoing wave at the event horizon. That is why the factor 2 appears in the power of the first 



multiplier in (3.61). The parameter tq is fixed to be the inner horizon r_, in order to move all the 



singularities outside the unit circle and, at the same time, to provide the best convergence of the 



infinite continued fraction (3.44). 



Since A can be calculated numerically as a function of uj (see section 3.4.5), we are able to solve 



the equation (3.44) with respect to u. 



In the Reissner-Nordstrom-de Sitter limit (a = 0) for neutral field (g = 0) we obtain K 



uor 



and 



R{r) 



r — 



r — r+ To 



(3.62) 



The tortoise coordinate is defined by dr^ 
event horizon. 



r^dr 



and K 



A' (r+) 

^ ^ — is the surface gravity on the 



3.4.7 Frobenius series for the massive scalar field equation in the 
higher-dimensional Reissner-Nordstrom-de Sitter background 

A D-dimensional Reissner-Nordstrom-de Sitter black hole is described by the metric 



ds"^ = f{r)df 



dr'^ 
W) 



r^dn 



D-2i 



(3.63) 
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where dQ£)-2 is the hne element of a (D — 2)-dimensional sphere, 

After separation of the angular and time variables, the radial part of the massive scalar field 
equation of motion (□ + IJ?)'^ = is reduced to the wave-like equation 

(3.65) 

where integer / = 0, 1, 2 . . . parameterises the angular separation constant. 
The Frobenius series for this case are 



i?(r)=(^ %-^(-V-^. r-^^^--^° ), (3.66) 

where k = -f'{r^), e*"^*-^-* describes the outgoing wave for the spatial infinity and the horizons 

^ j"^^ = ^ • The sign in the exponent is fixed by the quasi-normal boundary condition: 

or /(r) 

the real part of A{r Too) must be of the same sign as the real part of the eigenfrequency u. This 
choice of the sign makes the wave outgoing at the cosmological horizon (or spatial infinity). 
For massless field (/i = 0) this exponent is e*'^^* [20] and the Frobenius series read 

i?(r)= ('^V'e-r-^J^^:^^). (3.67) 
r -r-Q J \r - To Too - r+ ' 



Since for D = 4 we can choose tq = r_, we come to ( 3.62[ ) (s = 0). 



At the cosmological horizon we can observe the same asymptotical behavior of the exponent 



e'^^'''^ ~ e'^'-% r Too < oo. 



If A = 0, il* > 5, /(r) = 1 + o(r ^) we can write the Frobenius series as [18 



R{r) = ^ e-^\-^u ^ . (3.68) 
\r -tq J \r -tq J 

For D = 4: the term of order ~ r^^ in /(r) leads to the non-trivial contribution [57] 

iuj 

\r -To J \r -To J 

The same approach could be applied for the equations of motion for the Maxwell field and the 



gravitational perturbations, because their radial part can be reduced to the form (3.34) [5H1 159] . 
Practically, the continued fraction coefficients appear to be so complicated for such cases, that we 
are unable to compute quasi-normal modes with this method during reasonable time. 
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3.5 Horowitz-Hubeny method 



In order to find quasi-normal modes in the asymptotically anti-de Sitter space-times, usually we 
need to impose the Dirichlet boundary conditions at the spatial infinity. Thus, we can find the 
appropriate expansion for the function R{r) in (3.34) without consideration of the singularity point 

(3.70) 



at the infinity. This method was proposed by Horowitz and Hubeny [7]. Namely, we define 

r — r I 



where k is the surface gravity at the event horizon. 

If we substitute (3.70) into (3.34), we find that ^/'(z) satisfies the equation 

s{zW{z) + ^^'(z) + ""^m = 0, 

Z Z'^ 

where 



iVs iVt A^^ 

n=0 n=0 n=\ 

are polynomial functions of z. 

Since the function ipi^z^ is regular at the event horizon z = 0, we can expand it as 



anZ 



n=0 

The Dirichlet boundary condition at the spatial infinity V'(l) = reads 

oo 

= 0. 



n=0 



Substituting (|3.72|) into (|3.71|), we can find the recurrence relation for the coefficients a„ 



n— 1 



k=0 



n{n — l)so + ntQ 



(3.71) 



(3.72) 



(3.73) 



(3.74) 



The equation (3.74) allows to calculate all a„ starting from an arbitrary oq, which fixes the 
scale of ip{z). Substituting a„ into (3.73), we find the equation with respect to the eigenvalue u. 
Practically, since the series (3.73) are convergent, we sum over some finite number of terms and 
solve (3.73) with respect to u. In order to ensure that this truncation does not cause incorrect result 
we have to increase the number of terms until the value of u does not change within the required 
precision. 

Note, that, for the sum (3.73) to be convergent, r_ has to be a singular point of the equation 
(3.34), and all the other singularities, except r = r+ and r = r^o, must lay outside the unit circle 
|z| > 1. If both of these conditions are impossible to satisfy, we must use the continued fraction 
method with the appropriate fixing of the behavior of R at spatial infinity. 
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Chapter 4 



Perturbations of four-dimensional black 
holes 



4.1 Quasi- normal modes of Schwarzs child and 
Schwarzschild-de Sitter black holes 

Perturbations of four-dimensional black holes were studied extensively in the context of possibility 
to observe quasi- normal ringing with the help of gravitational wave detectors (see reviews [SU [60]). 
Since in our Universe the cosmological constant appears to be positive, its correction to the quasi- 
normal spectrum must be taken into account. Gravitational perturbations of spherically symmetric 
black holes in the de Sitter background were studied for the first time by the numerical integration 
(see sec. 



2.2) in pO]. 



In |9] we study for the first time quasi-normal modes of massless Dirac and electromagnetic 
field in the Schwarzschild-de Sitter background. Also we study massless scalar and gravitational 
perturbations of Schwarzschild-de Sitter black holes, using the 6-th order WKB approach (see sec. 



3.3), which provides more accurate results. Namely, we substitute the effective potentials for the 



scalar, electromagnetic, Dirac and gravitational perturbations (2.11, 2.13, 2.14, 2.16, 2.17) into the 



WKB formula (3.33) and find, that the presence of the cosmological constant leads to decrease of 



the real oscillation frequency and to a slower decay (see fig. 4.1 ). 
For large / the following analytical expressions were found [H] 

VI - 9M2A 



Vl - 9M2A 
3v^M 



n + 



+ 



+ 



1 



(4.1a) 
(4.1b) 



As the black hole mass approaches its extremal value, the effective potentials (2.11, 2.13, 2.14, 2.16 



2.17) look like the Poshl- Teller potential (see sec. 3.2). Thus, using the Poshl- Teller approach, 
it was found in |15] that in the near-extremal limit of the Schwarzschild-de Sitter black hole the 
quasi-normal frequencies are proportional to the surface gravity at the event horizon k = ^/'(r+) 




(4.2a) 



(4.2b) 
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Figure 4.1: Real and imaginary parts of the fundamental quasi-normal frequency of the 
Schwarzschild-de Sitter black hole for scalar (Z = 2), electromagnetic (/ = 2), gravitational (/ = 2), 
and Dirac (k = 2) perturbations plotted with solid blue, dashed red, dot-dashed green and dotted 
black lines respectively. 

Despite the dominant contribution of the imaginary part does not depend on the field spin s, 
the metric perturbations decay slower among the considered perturbations for all values of A. 



4.2 High overtones in the quasi- normal spectrum 

Though, practically, we can observe only low-damping oscillations, the frequencies with large imag- 
inary part attracted some attention in the context of the black hole thermodynamics. Such frequen- 
cies are called high overtones. Since their imaginary part is very large, the WKB approach is not 



appropriate and one should use the more accurate Frobenius method (see sec. 3.4). 

It has been suggested in [61] that the asymptotic value of the real part of the quasi-normal 
frequency (i. e. when the imaginary part approaches infinity) coincides with the so-called Barbero- 
Immirzi parameter [U^ This parameter must be fixed in order to predict the Bekenstein- 

Hawking formula for entropy within the framework of Loop Quantum Gravity. For Schwarzschild 
black holes it was found numerically [5l] and analytically [Mj that the asymptotical behavior of the 
quasi-normal frequency is given by 

for the gravitational perturbations and for the test scalar field, while for the electromagnetic per- 
turbations the real part asymptotically approaches zero 

We have shown numerically in [TT] that for perturbations of the massless Dirac field, the asymp- 
totical value of the real part of the highly damping quasi-normal frequencies is also zero (see fig. 



4.2). The correct spacing of the imaginary part was found later in |65], with an alternative form of 
the effective potential for the Dirac perturbations. The asymptotical formula for the high overtones 
of the massless Dirac field in the Schwarzschild background reads 



m 

^nr+ — . 
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Figure 4.2: Real part of the highly damping quasi-normal frequencies as a function of imaginary 
part {k± = 1 perturbations of the massless Dirac field in the Schwarzschild background). The 
frequencies are measured in inverse units of the black hole mass. 



Using the Frobenius method (see sec. 3.4), for Schwarzschild-de Sitter black holes we have found 
for the first time in pU] that again the real part of the quasi-normal frequencies for electromagnetic 
perturbations asymptotically approaches zero, satisfying 



iKn 



U 



LOr, 



where k = \f'{ 



and hie 



are the surface gravities at the black hole horizon and the 



cosmological horizon respectively (see fig. 4.3). This result was later confirmed analytically in |66j . 

The real part of the quasi-normal modes for metric perturbations does not approach a constant. 
Frequencies with high imaginary part satisfies the non-algebraic equation 



cosh 



-nuo nuj 



3 cosh I 1 



0, 



which implies oscillation of the real part as a function of the imaginary part (see fig. 4.4). 



4.3 Decay of charged scalar and Dirac fields in the Kerr- 
Newman-de Sitter background 

For charged black holes, the scalar field electrodynamics can describe the interaction of a charged 
field with the electromagnetic background of the black hole. When the influence of the spin of 
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Figure 4.3: Real part of highly damping quasi-normal frequencies as a function of imaginary part 
for Schwarzschild-de Sitter black hole (/ = 1 electromagnetic perturbations, AM^ = 0.02). The 
frequencies are measured in inverse units of the black hole mass. 



the field is neglected, we can consider the charged scalar field. The late-time tails of the charged 
scalar hair were studied for the first time by Hod [671168] in the context of the gravitational collapse. 
Quasi-normal modes of the charged scalar field for various black holes were studied for the first time 
in [69j. Quasi-normal modes of the charged massive scalar field for Reissner-Nordstrom black holes 
were considered in [7D]. Decay of the charged Dirac field in the Reissner-Nordstrom background 
was studied in [7T]- The decay law for the late-time tails of the charged massive Dirac field also was 
found for Reissner Nordstrom [72] and Kerr-Newman [73] black holes. Yet, the calculations of the 
quasi-normal modes in the above papers were limited by the third order WKB accuracy. The WKB 
accuracy does not allow to study the quasi-normal modes of the near extremal charged black holes. 
Namely, the decay rate of charged fields in the near extremal Reissner-Nordstrom background is 
the same within the WKB accuracy. 

The numerical study of quasi-normal frequencies of the charged fields (both scalar and Dirac) 
in the Kerr-Newman-de Sitter background was done in [22] within a much more accurate Frobenius 



method (see sec. 3.4.5, 3.4.6). We have shown that for not very large value of Q, the charged field 



decays quicker than the neutral one. For the near extremal value of Q, the charged field decays 



slower than the neutral one (see fig. 4.5). 



Let us summarise briefly what happens if a black hole has charge and rotation. We measure all 
the quantities in units of the size of the event horizon. 

• As we know from previous sections, the presence of the cosmological constant decreases the 
absolute values of the real and imaginary parts of uo. 
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Figure 4.4: Real part of highly damping quasi-normal frequencies as a function of imaginary part 
{1 = 2 metric perturbations of the Schwarzschild-de Sitter black hole, AM^ = 0.02). The frequencies 
are measured in inverse units of the black hole mass. 
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Figure 4.5: Real part and imaginary part of the fundamental quasi-normal frequency of charged (g) 
scalar field (/ = 0) for the charged (Q) black hole. 
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We already know, that as the multipole number I increases, the real part of uj grows, while 
the imaginary part of u approaches a constant. 



Metric perturbations decay slower than any field perturbations of the same /. 



The dependance of u on charge of the field q and of the black hole Q is shown on figure 4.5 



— If qQ > 0, the real part of uj monotonically grows as the field charge q increases. It grows 
also as the black hole charge Q increases, attains some maximum value at a large (close 
to extremal) value of Q, and then decreases. 

— If qQ < 0, the real part of u decreases when either the black hole charge or the field 
charge increases. 

— The imaginary part of uj has more difficult behavior, in most cases decreasing its absolute 
value as qQ grows. 

The influence of the rotation parameter depends significantly on the projection of the field 
momentum on the axis of the black hole rotation m. We will consider this dependance later 



for the massive scalar field in the Kerr background (see sec. 6.3). Qualitatively the same 



behavior was observed for Kerr black holes, projected on the brane (see sec. 5.3.1). 
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Chapter 5 



Perturbations of higher-dimensional 
black holes 



5.1 Stability and quasi- normal modes of Reissner- 
Nordstrom-de Sitter black holes 



Last years higher-dimensional black holes have attracted considerable interest in the context of 
the string theory and higher-dimensional brane- world models [TH [7S1 [THl ETJ [7H1 U^- Some of 
such models allow compactification radius of extra dimensions to be of macroscopic size [501 1ST] . 
That is why the extra dimensions are important if we study small black holes. Such black holes 
could be produced at the next-generation particle colliders, probably at energies of order ~ ITeV 
[821 [831 Ha Eg. 

Being an important characteristic of a black hole, its quasi-normal spectrum can be used in 
future experiments to find number of extra dimensions and other parameters of the theory. That is 
why we study quasi-normal spectrum of black holes within different scenarios with various possible 
numbers of extra dimensions. 

When a black hole is much smaller than the size of extra dimensions, it can be described by 

dimensional Einstein theory. The D- dimensional generalisation of the Schwarzschild metric was 
done by Tangherlini [5S] . 

Let us consider the D-dimensional Reissner-Nordstrom-de Sitter black hole, which is described 
by the line element 

fir)dt'-^-r'dnl (5.1) 



where dfld is the line element of ci-sphere, d 
quantities X, Y, Z are defined as 



X 



dMAa 



Y 



fir) 
D-2, f{r 

2Ar2 
d{d+l)' 



1 — X + Z — Y . The dimensionless 



r2d-2 ■ 



where Ad 



27r('i+i)/2 



r((d + i)/2) 

hole charge and A is the cosmological constant. 



is the area of a unit (i-sphere, M is the black hole mass, Q is the black 



Perturbations of the Einstein-Maxwell equations can be reduced to the wave-like form (2.10). 



The corresponding effective potential depends on their transformation law under rotations on the 



o?-sphere (2.15). Thus, there are three types of perturbations described by three effective potentials: 
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tensor (Vt) 

EH] 



vector (VV) and scalar (V^) type. The explicit form of the potentials was derived in 



Vrir) 
Vvir) 



fir) A , ^ , f'{r)rd ^ f{r)d{d-2) 



(5.2) 



\+d±\ X 



+ 2ZXd{d-l)x 



d{d - 2)(1 -Y) + Zd{5d - 2) - X{d^ + 2) 



Vs{r) = fir) 



(5.3) 
(5.4) 



where 



1 - 



A + 



d{d+l\ 
2 

d{d+l 



■5X, 



(1 + X5)X, 



[- 



-Ad^d + 2){d+ 1)^6^X^ - A8d\d + l){d- 2)6X - 16{d - 2)(rf - 4)] Y 
S^d\3d -2){d+ l)\l + X6)X^ + 46^d\d + 1)2 {{d + l){3d - 2)X6 + Ad^ + d-2}X 
+A5{d +l){{d-2){d- A){d + 1)(A + d^)6 - 7d^ + 7d^ - lAd + 8} 

2): 

2)A(1 + X6)X - 16(rf -2){d- A)X^] Y 



-Ad^d 

J3/ 



+ {I6{d + 1) (-4A + 3d^{d - 2)) 5 - 16(3rf - 2)(rf - 2)} X + 64A + 16d{d + 2) 

^ .,'^/, „\/, ^N9/^ .,,-\9-r^9 . „ ,0 , , , , 



2){d+l)\l + X6)^X^ + 48d\d+l){d , ...... 

d^{3d -2){d+ 1)^5(1 + X6)^X^ - Ad\d + 1)2(1 + A5)2 {{d + l){3d - 2)X6 - d^} X^ 
+A{d + 1)(1 + X6) {X{d -2){d- A){d + 1)(A + d^)6 + Ad{2d^ -3d + 4)A 
+d'^{d -2){d- A){d + 1)} + 64A^ + 16d{d + 2) A' 



16A {(rf + 1)A (-4A + 3d'^{d -2)) 6 + 3d{d - 4)A + 3d\d + l){d - 2)} X . 



Here we defined X = {I + d){l — 1) through the multipole number / = 2,3,4 
dimensionless constant 



The value 5 is a 



2A5 



1 + 



16AZ 



(rf+ 1)2X2 



1. 



Note, that the effective potential for tensor-type perturbations coincides with the potential for 



a test scalar field (3.65) for / > 2. Vector-type and scalar-type metric perturbations are coupled to 
perturbations of the electromagnetic field. For neutral {Q = 0) black holes the type potentials 
are reduced to the effective potential for the test Maxwell field and for the type "-" Vy and 
describe pure vector-type and scalar type metric perturbations respectively. 

It was proven analytically that the black hole is stable against tensor-type and vector-type 
gravitational perturbations |58]. The effective potential of the scalar type (5.4) is not positive 



definite, having a negative gap for small /. Since the effective potential has an extremely complicated 
form, the appropriate ansatz for the S-deformation technique (see sec. 2.3.2) was not found. That 



is why, the black hole stability against scalar-type perturbations was not proven analytically. 

The stability of higher dimensional Schwarzschild-de Sitter black holes was proven numerically 
in [20]. We have shown in time domain that the gravitational perturbations of scalar type decay 
for arbitrary black hole mass and A-term (see figs. [5lp:2| . 
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Figure 5.1: Effective potentials and time- domain profiles for scalar-type gravitational perturbations, 
D = 5 (blue). . . D = 11 (red) (Z = 2, Q = 0, A = 0). For higher D both the peak and the negative 
gap of the potential increase. Profile for higher D decays quicker. All quantities are measured in 
units of the event horizon r+. 
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Figure 5.2: Time-domain profiles for scalar-type gravitational perturbations of 11-dimensional 
Schwarzschild-de Sitter black hole {Q = 0, D = 11, p = r+/roo) for p = 0.3 (blue), p = 0.5 
(green), p = 0.7 (yellow), p = 0.8 (orange), p = 0.9 (red). Profile for higher p decays slower. All 
quantities are measured in units of the event horizon r_|_. 
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Figure 5.3: Time-domain profile for gravitational perturbations of scalar type "-" of the Reissner- 
Nordstrom-de Sitter black hole {D = 11, p = r^/voo = 0.8, / = 2) for various values of the black 
hole charge q = Q/Qext- q=0.4 (brown) q=0.5 (blue) q=0.6 (green) q=0.7 (orange) q=0.8 (red) 
q=0.9 (magenta). The smaller q, the slower growth of the profile is. 



Also we studied the quasi-normal modes for all types of gravitational perturbations of D- 
dimensional spherically symmetric black holes. Since for the observational purposes only the dom- 
inant quasi-normal frequency is essential we were limited by the fundamental mode (/ = 2) only. 
However, using the WKB approximation, we are able to find the large / formula for any given black 
hole parameters. 

The dependance of the quasi-normal modes on the A-term is qualitatively the same as for the 
four- dimensional black holes: both real and imaginary parts of the quasi-normal frequency decrease 
their absolute values as the cosmological constant grows (see fig. 5.2). We also observe that the 
real and imaginary parts of the quasi-normal frequency are non-monotonic functions of the black 
hole charge, in most cases decreasing their absolute values as Q grows. 

In we have shown for the first time that Reissner-Nordstrom-de Sitter black holes are 
gravitationally unstable for large values of the electric charge and cosmological constant in D > 7 
space-time dimensions. On the figure 5^ we see the time-domain profiles for the linear gravitational 
perturbations of scalar type "-" (5.4) of the near extremal Reissner-Nordstrom-de Sitter black hole 
in D = 11 dimensions. For sufficiently small values of the black hole charge we observe usual 
picture of the quasi-normal ringing. Then, as the black hole charge increases, a purely imaginary 
(non-oscillating) mode becomes dominant, decreasing its decay rate until the threshold point of 
instability is reached. After crossing the instability point we observe the growing non-oscillating 



mode (see sec. 2.3.2). Its growth rate increases as the black hole charge grows. Therefore, we 
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Figure 5.4: The parametric region of instability in the right upper corner of the square (p = r+/roo, 
q = Q/Qext) for D = 7 (top, black), D = 8 (blue), D = 9 (green), D = 10 (red), D = 11 (bottom, 
magenta). 

conclude that exactly at the threshold point of instability there is some static solution {u = 0) of 
the perturbation equation. The static solution at the threshold point of instability was observed 



also in time-domain for the black strings (see sec. 6.4). 



The parametric region of instability is shown on the figure 5.4 The larger number of space 



time dimensions D is, the bigger region of instability we observe. Though the region of instability 
increases with D, the charged black holes in the asymptotically fiat space-time are stable at least 
for D < 11. The instability occurs if both the black hole charge and the cosmological constant are 
large enough. 



5.2 (In) stability of L)- dimensional black holes in the Gauss- 
Bonnet theory 

Higher dimensional quantum gravity implies corrections to classical general relativity. The dominant 
order correction to the Lagrangian is called the Gauss-Bonnet term [87j. This term is squared in 
curvature and vanishes for D = 4. The effective action is given by 

S = j d'^x^g {R + aiRabcdR'''''' - ^RcdR"" + R^)) , (5.5) 

where a is a positive coupling constant, [a] = L^. 
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The spherically symmetric black hole solution, which satisfies the corresponding equations of 
motion, is described by the line element (5.1) with 



,2 



which is reduced to the Tangherlini metric [HE] in the limit of a ^ 0. 

Quasi-normal modes and late-time tails of the test scalar field for black holes in the Gauss-Bonnet 
theory were studied in [881 ES] . 

The quantum corrections imply that the effective potential for tensor-type gravitational per- 
turbations (Vt) does not coincide with the potential for the test scalar field [90]. The effective 
potentials for vector-type (VV) and scalar-type (V5) gravitational perturbations were found in 
They are given by the formulae 



(D-2)(c+l) / B{r)\ 1 d- 



VT{r) = firy- 3 + , ^ — :j-Wr''-'Air)qir), (5.7) 



r2 V Mr) J v^r^-M(r)g(r) dr^ 



f{r)U{r) 

64r2(L) - 3)2A(r)2g(r)8(4cg(r) + {D - l)R{q{rY - 1))^ 



where we used the following dimensionless quantities 



A{r) 
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g(r)2 \2 D-3J \2 D-3 



B{r) = A{rf[l + -^] + (l ^ 
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D-2 ^' ^ = 2'^'^ 
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f/(r) = 5{D-lfR\l + R)-3{D-lfR{{D-l)R'^ + 24c{l + R))q{r) + 

+2{D - l)^(24c(D - 1)R^ + 168c2(l + R) - {D - l)R\-3 + 5R + 7D{1 + R)))q{rf + 

+2{D - l)^i?(-184c^ + iD- 1)(13 + D)R^ + c(-84 + 44i? + 84D(1 + R)))q{rf + 

+ {D - l)3(384c^ - 48c(2 + D{3D - 5))R^ + 192c2(-ll + D + (-15 + D)R) + 

+ {D - l)i?2(-3(7 + 55i?) + D{2Q + lOQR + 7D{1 + R))))q{r)^ + 

+ {D - lfR{-64:C^{D - 38) + (L> - 1)(71 + D{7D - 90))R^ + 

+16c(303 + 255i? + 13D^{1 + R) - 2D{73 + 81R)))q{rf + 

+4{D - l)2(96c^(-7 + D)- 8c{D - 1)(145 - 74D + QD'^)R^ - 

-8c^(9 - 175i? + D(-58 - 34i? + 11D{1 + R))) + {D - l)R^{-5{79 + 23i?) + 

+^(5(57 + 41i?) + D{-81 - 89R + 7D{1 + R)))))q{rf - 

-A{D - iyR{8c^{43 + (72 - 13D)D) + {D - l)(-63 + D(99 + D(-49 + 5D)))R^ + 
+4c(321 + 465i? + D(121 - 39i? + D{-123 - 107R + 17D{1 + i?)))))g(r)^ + 
+ {D - l)(128c3(-9 + D){D - 5) + 32c(D - 1)(246 + D{9 + L>(-55 + 8D)))R^ + 
+64c\D -5){D^-3 + (49 + {D - A)D)R) - 

-{D- l)i?2(ii73 + 565i? + D(-4(997 + 349i?) + ^(6(393 + 217i?) + 

+D(-548 - 452i? + 45L>(1 + R))))))q{rf + 

+ {D - l)R{-6Ac^{D - 5)(36 + L>(-13 + 3D)) + 

+ {D - 1)(635 + D(-1204 + 3L>(294 + D{-92 + 9D))))R^ - 

-8c{D - 5)(63 + 31i? + L>(127 + 191i? + L>(-47 + D + (-79 + D)R))))q{r)^ + 

+2{D - 5)(64c^(D -5){D-3) + 8c{D - l)(-27 + D(141 + (-43 + D)D))R^ + 

+8c^{D - 5)(-3 + 77i? + D{D - 2 + {D - 18)R)) + {D - 1)^R'^{-33{R - 7) + 

+^(59 + 43i? + L>(-59 - 35i? + 9L>(1 + R)))))q{r)^° - 

-2{D - 5)i?(24c2(-ll + D){D - 5){D - 3) + {D - 1)\-Q5 + D{81 + D{7D - 39)))R^ + 

+12c(-7 + D){D- 5)iD -3){D- 1)(1 + i?))g(r)" + 

+ {D - 5)^(-l + D)R\lQc{2Q + {D - 9)D) + 

+ (D - 1)(77 -3R + D{-18 + D + {D- 2)R)))q{r)'^ + 

+ {D - 5)\D - 3)\D - l)2i?3g(r)i3. 



The quasi-normal modes of gravitational perturbations of Einstein-Gauss-Bonnet black holes 



were found in time-domain (see sec. 3.1) for the first time in |23]. Consequently, we confirm the 
instability of these black holes in five- and six- dimensional space-times, proven analytically in |92j . 

Namely, in five dimensional case scalar-type perturbations are unstable for a > 0.207r^ while 
in six dimensions tensor-type perturbations are unstable for a > l.OOGrl. We can see on the figure 



5.5 that the larger /, at the earlier times instability growth occurs, and the stronger the growth rate 
is. In the region near the threshold value of a, one can observe the growth only for large enough I, 
while perturbations of lower multipole number are not growing. For each / there is some maximal 
value of a for which the perturbations are not growing. In order to find the threshold value of a 



numerically we extrapolate this value for / — * oo (see fig. 5.6). Einstein-Gauss-Bonnet black holes 
in D > 7 space-times are stable. 

In units of the event horizon radius, the imaginary part of fundamental quasi-normal modes de- 
creases, when a increases, for all numbers of D and all types of perturbations. Unlike the imaginary 
part, the real oscillation frequency does not behave uniformly: it decreases as a grows for most cases 
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Figure 5.5: The picture of instability of tensor-type of gravitational perturbations of Gauss-Bonnet 
black holes, developing at large multipole numbers: D = 6, I = 8 (red), I = 12 (green), / = 16 
(blue), a = 1.3. All quantities are measured in units of the event horizon r+. 

of tensor and vector modes. The behavior of the scalar mode is different: there are two competing 
for the domination modes at different stages of the quasi-normal ringing. This superposition of 
modes, also with competing excitation coefficients, makes dependence of the fundamental scalar 
type quasi- normal modes on a and D non-monotonic. On the figure [STT] one can see that at the first 
stage the actual frequency of the dominant mode is much larger than at the second stage, while 
their damping rate stays almost the same. 



5.3 Quasi- normal modes of brane- localised Standard 
Model fields 

The other possible scenario is that the Standard Model particles (scalars, fermions and gauge bosons) 
are restricted to live on the 3 + 1-brane, which is embedded in the higher-dimensional bulk, while 
the gravitons can propagate also in the bulk. Therefore, if we study propagation of the fields near a 
D-dimensional black hole, we must consider the induced-on-the-brane gravitational background. If 
the size of extra dimensions is large, comparing to the size of the black hole, the induced background 
is given by the projection of the £)-dimensional black hole metric onto the brane by fixing the values 
of the additional angular coordinates that describe the {D — 4) extra spacelike dimensions [931 El] • 
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Figure 5.6: The threshold a as a function of the inverse multipole number / for tensor type of 
gravitational perturbations of Gauss-Bonnet black holes D = 6. The points I = 16, 20, 32, 40, 50, 64 
were fit by the line a = 2.627/"^ + 1.005. The theoretical result is at ~ 1.006 (the value of a is 
measured in units of the event horizon r+). 



5.3.1 Kerr black holes 

The quasi-normal ringing of brane-localised fields propagating in Schwarzschild, Reissner- Nordstrom 
and Schwarzschild- (anti) de Sitter induced gravitational backgrounds was studied in [95]. Yet, 
during the high-energy collisions of elementary particles resulting in the creation of black holes, it 
is unnatural to expect that only head-on collisions, leading to spherically symmetric black holes, 
would take place. Collisions with a non-vanishing impact parameter are most likely to occur, and, 
in addition, it is for these collisions that the black- hole production cross-section is maximised [96]- 
Therefore, microscopic rotating black holes should be the most generic situation, and the effect of 
the angular momentum of the black hole on the quasi-normal spectra of brane-localised fields is 
essential and cannot be neglected. This effect has been studied in [T7]. 

The line-element, describing a higher-dimensional rotating neutral black hole, is given by the 
Myers-Perry solution [UT]. After the projection of the Myers-Perry metric onto the brane, the brane 
background assumes the form [HS] 

= y--^;^)'^^ + dtdip--dr -lldO -[r +a + ^^^^^ J sm 9 , (5.10) 



where 

A = -I- 



r^ + a^-^, J: = r' + a'cos'9. (5.11) 



The parameters /i and a are related to the mass M and the angular momentum J, respectively, of 
the black hole through the definitions 



(D - 2) 7r(^-i)/2 2Ma 
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Figure 5.7: The picture of time-domain evolution for scalar-type gravitational perturbations of 
Gauss-Bonnet black holes D = 10, I = 2, a = 0.01. One can see that two modes are dominating at 
the different stages. All quantities are measured in units of the event horizon r+. 



where = SirG = Sir/M^ ^ is the D-dimensional Newton's constant. The radius of the event 
horizon parameterises the black hole mass 



fx 



We should note here that the higher-dimensional black hole is assumed to have only one non- 
vanishing component of angular momentum, about an axis in the brane. This is due to the sim- 
plifying assumption that the particles that created the black hole were restricted to live on an 
infinitely-thin brane, therefore, during collision they had a non-vanishing impact parameter only on 
a 2-dimensional plane along our brane. 



The equations of motion of the Standard Model fields in the background (5.10) can be reduced 
to the following to coupled equations for the angular 
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and radial parts 
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Figure 5.8: Fundamental quasi-normal modes for the 6-dimensional black hole projected on the 
4-brane. 



where A is the separation constant, m is the azimuthal number, s = 0, 1/2, 1 for scalar, fermion and 
gauge boson fields respectively and K is defined as 

K 



+ a") 



am. 



In order to solve numerically the equations (5.13) and (5.14) we use the Frobenius method. Since 



(3.55) coincides with (5.13) in the limit of a — 0, the equation for the angular part can be solved 



as described in the section 3.4.5 The analysis of the singular points of the equation (5.14) allows 



to find the Frobenius series for the function R{r) [T7] 



R{r) 



r — Ta 



r — tq \ r — To 



-1/3 oo 
i=0 



r — Ta 



r — Tq 



(5.15) 



where (3 is fixed (see sec. 3.4.6 ) as 



ti;r+(r^ + a^) — mar^ 



K{r+) 

A'(r+) (L)-5)(r2 +a2) + 2r2 



From the figure |5^ one easily observes that, as a increases, the absolute values of both the real 
and imaginary parts of the fundamental quasi-normal frequency decrease. This makes the damping 
time longer and the field oscillations on the brane longer-lived. For large a one can observe that the 
lifetime of a fermion signal is longer for the positive value of the azimuthal number m. The same 
behavior was observed also for the scalar field: the field oscillation lifetime grows as m increases. 



On the figure 5.9 , we display a few of the higher overtones for a brane-localised gauge field (with 
/ = 3, m = 0). For all higher overtones, it was found that an increase in a leads again to the 
decrease of the absolute values of both the real and imaginary parts. 
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Figure 5.9: A few higher overtones of the electromagnetic field for the 6-dimensional Kerr black hole 
projected on the 4-brane (/ = 3, m = 0) for the range (0,r_|_) of the angular momentum parameter 
a with a step of r+/8. 

For fixed a we find that, as D increases, the imaginary part of the quasi- normal frequency 



decreases while the real part changes insignificantly (see fig. 5.10). This behavior of the real part is 
different from the case of dimensional black holes, which are not projected on the brane (see sec. 



5.1). For those black holes the real part of the quasi-normal frequency also grows as D increases. 



Therefore, we conclude that higher dimensional black holes are better oscillators than the black 
holes, projected on the 4-brane. 



5.3.2 Gauss-Bonnet black holes 

Within large extra dimensions scenarios of TeV-scale gravity, the classical space-time, induced by 
mini black holes, has large curvature along the transverse collision plane. Thus quantum gravity 
effects, and in particular higher curvature corrections to the Einstein gravity, cannot be ignored 
[98j . Therefore, the quantum corrections, provided by the Gauss-Bonnet theory (see sec. 5.2), must 
be taken into account, when mini black holes are considered. 

In this context, the propagation of brane-localised Standard Model fields in the background, 
induced by the Gauss-Bonnet black hole, has been studied numerically for the first time in [25] . 

The metric induced on the brane is given by the following line element 



(5.16) 



with the function /(r) defined in (5.6). 
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Figure 5.10: Fundamental quasi-normal modes for the higher-dimensional Kerr black hole (a = r+) 
projected on the 4-brane. 



The effective potentials for the Standard Model fields can be derived in the same way as in the 
Schwarzschild-de Sitter case. Their explicit form can be obtained by substituting the function /(r) 



from (5.6) into the effective potentials for the test scalar (2.11), massless Dirac (2.13), and Maxwell 



(2.14) fields. 



As well as for the Schwarzschild black holes, the oscillations of brane-localised Standard Model 
fields decay faster for higher D. On the other hand, the Gauss-Bonnet term causes the perturbations 
decay slower. The real part of the quasi-normal frequencies has a more complicated behavior: for 
D = 5 it decreases as a grows, but for higher-dimensional cases it starts growing first and then 



decreases after some value of a is reached. However, the quality factor Q 



grows for all fields and all values of D (see fig. 5.11). 



Re(uj) 
2Im((^) 



mcreases as a 



As we can see in time-domain (see fig. 5.12), the late-time tails decay according to the inverse 
power law, which is found to be [23] 

^ocr(2'+^-i). (5.17) 

This law depends only on the multipole number / and the number of extra dimensions D. 

Using the WKB formula, we are able to find the Gauss-Bonnet corrections to the large multipole 
limit 



uj = fiij(l + aAi + . . 

LJ = n^il + aAi + . . - i^i (1 + aBi 



O 



n 
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Rq 
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Figure 5.11: The quality factor of the scalar field localised on the 4-brane as a function of the 
Gauss-Bonnet parameter a {D = 9, I = 1). 

where Rq is the point, where the effective potential reaches its maximum for large / 



Rq 



D-3 



D - 1 



O 



We see that for large multipoles the quality factor also decreases as D grows 

O ^ 

The corrections of the first order of a are given by 

1 D - 4 fD - 1 r 



^1 



r\D-\\ 2 
\ D - 4 ( D -\ 



rlD-1 



{D-2 



Rn 



< 0, 



implying quicker decreasing of the absolute value of the imaginary part than decreasing of the real 
part. It means that the quality factor grows as a increases. Thus we conclude, that the large 
multipole limit resembles the main properties of the fundamental quasi-normal modes for small 
multipole number. 
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5 10 50 100 500 1000 

Figure 5.12: The time-domain profiles of tlie brane-localised massless scalar field for the Gauss- 
Bonnet black hole (D = 7, a = 5) for / = (blue), / = 1 (green), / = 2 (orange), I = 4 (red). The 
bigger / corresponds to the longer life of quasi-normal ringing and the quicker tail decay. 



5.4 Perturbations of squashed Kaluza-Klein black holes 



If the compactification radius of the extra dimensions is comparable with the size of the black 
holes we must take into account the size of extra dimensions. Such higher dimensional model of 
black holes has the asymptotic structure of the Kaluza-Klein type. The simplest example of five- 
dimensional black objects with the Kaluza-Klein geometry is the black string, the direct product of 
four-dimensional black hole and a circle. These objects look different from four-dimensional black 
holes only at sufficiently high energies, when Kaluza-Klein modes are excited. Therefore within these 
space-times we need high energy regime to see the extra dimensions. At the same time, there exist 
exact solutions of Kaluza-Klein black holes with squashed horizons, that look like five- dimensional 
squashed black holes near the event horizon, and like a Kaluza-Klein space-time at spatial infinity. 
Owing to the non-trivial bundle structure, the size of the extra dimension might be observed even 
at low energies by detecting e. g. their Hawking radiation |99j . 

To the best of my knowledge, quasi-normal frequencies of such black holes have been studied for 
the first time in [M] . 
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5.4.1 Quasi- normal modes of the scalar field for rotating squashed 
Kaluza-Klein black holes 

The five-dimensional rotating squashed Kaluza-Klein black hole with two equal angular momenta 
is described by 



ds^ = de - ^k{rfdr^ 

^0 



[k{r){dal + dal) + dal 



{dt--da^f, (5.19) 



+ d 



with 



dai = — sin ipd6 + cos i/j sin 6d(j) , d(72 = cos il)d6 + sin if) sin 6d(j) , da^ = dip + cos 6d(j) , (5.20) 
where 0<^<7r, O<0<27r and < ip < Att. The parameters are given by 



So 
Ao 

k{r) 



(r^ + a^f - Mr"^, 



(5.21) 



Here M and a correspond to mass and angular momenta, respectively. Values r = r+ and r = r_ 
are outer and inner horizons of the black hole. They relate to M and a by = {r^r^y,M — 
2a? = r"^ + r"^. The parameter r^o corresponds to the spatial infinity. In the parameter space 
< r_ < r+ < r^o, r is restricted within the range < r < Too- The shape of black hole horizon is 
deformed by the parameter k{r^). 



by 



The wave equation for the massless scalar field ^{t,r,6,(f),ilj) in the background (5.19) is given 



(5.22) 



Taking the ansatz 



-iuit 



R{p)e 



imil>+i\il} 



S{9), 



where >S'(^) is the so-called spheroidal harmonics, the radial and time variables can be decoupled 
from angular ones, so that the final wave-like equation reads 



d 

dp 



dR{p) 
dp 



+ 



+ A-/(/ + 1) + A2 



R{p) = 0, 



where / is the non-negative integer multipole number, \m\ < I and |2A| < 21 are integers. 
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Figure 5.13: Real and imaginary part of the fundamental quasi-normal frequency of the test scalar 
field (/ = 0) for rotating squashed Kaluza-Klein black holes. 



The radial coordinate p is given by 



P = Po-2 ^, (5.27) 



r^-r 



with 



2 ^0 / 2 I 2\ 

PO = -ji^oo+a ), 

. ,(, = 0) = ^^^±^^^^^. (5.2S 



Note that the three parameters po and p± = po^±/(^^ ~ ''^i) '^^^ define the metric (5.19) if 
Too < oo. In some papers they are used to parameterise the black hole, instead of the parameters 
r-oo, r±. 

The quasi-normal modes have been found with Frobenius method (see sec. 3.4), using the 
following expansion [21] 

9 9 \ —il3 oo / 9 9 \ n 



''={7^} ^""p-"-'Y.-''[7^) . (6-29) 

n=0 



where (3, v and Vt are chosen in order to eliminate the singularities at r = r+ and r = Too. The sign 
of a and f2 is chosen in order to remain them in the same complex quadrant as uo. 



The fundamental quasi-normal modes are presented on the figure |5.13[ One can see that the 
real oscillation frequency exerts some irregular growth (with local minimums) when Too is increasing 
until some moderately large values of Too. At larger Too the growth of Re(co') changes into monotonic 
decay. The imaginary part of uj that determines the damping rate also has some initial irregular 
growth when increases, but at larger r^o the two scenarios are possible: either monotonic decay 
(for large values of r_) or monotonic growth (for small and moderate r_). Thus, for a given mass 
and angular momentum of the black hole, one can learn the size of extra dimension Too from values 
of quasi-normal modes of the emitted radiation. 
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5.4.2 Gravitational quasi-normal modes for non-rotating squashed 
Kaluza-Klein black holes 

In the previous subsection, we have considered the scalar field in the background of squashed Kaluza- 
Klein black holes. However, the tensor perturbations are more interesting from the point of view of 
the stability and possibility to observe gravitational waves from black holes. 

The metric of the uncharged non-rotating squashed Kaluza-Klein black hole is a particular case 
of the metric (5.19) 



ds^ = Fip)dT^ - ^^rfp' - Ap'Gipfda+da 



3\2 



;{da') 



(5.30) 



where we have defined r = 2pQt/roo and 



F{p) = 1 
Here, we have used a basis 



G{p)' 



1 + 



Po 



4po(p+ + Po) 



(icT^ = -{da^ =F ida'^) . 



Since the space-time has the symmetry SU{2) x f/(l), the metric perturbations can be classified 
by eigenvalues J, M for SU{2) and K for U{1). Here we consider only zero modes J = M = 0. Even 
in this case, since da^ carry eigenvalues K = ±1, each component could have different eigenvalue 
K. It is important to recognise that the components with different K are decoupled. That is why we 
have the master equations for each K. To obtain master equations, we choose the gauge condition 
as 



h 



t3 



0. 



(5.31) 



As is shown in |lUOj . the perturbation equations for the = 2 mode can be reduced to the wave 
equation for with the effective potential in the form 



-{p+- p) 



64p5 + 256pVo - 32p3po(p+ - Ilpo) + ^P+pl{p+ + Po) 



16p3po(p+ + Po)(p + Po)^ 

+8pVo(2p+ - 5p+po + 25p^) + pPo(20p^ - 9p+po + 35p2) 



(5.32) 



Similarly, the perturbation equations for the \K\ 
for /ip+ with the effective potential 



1 mode can be reduced to the wave equation 



P+\ \^ ^ 7(-p+ + p) ^ 17p+ - 19p ^ 9(-3p+ + 7p) ^ P+-P 
p / [p+po 16(p + po)^ 8p(p + po)2 16p2(p + po) p\p + p+ppq 
Kp+-Pfp 2(p+ + 2p) 



(P+ - 2p)(p+Po - p(p + 2po))2 (p+ - 2p)(p+po - p(p + 2po)) 
Finally, for the K = Q mode, we obtain the wave equation for /133 with the effective potential 

-(P+ -P) 



(5.33) 



16p3(p + po)3(4p + 3po)2 



256p+p' + 64p^(17p+ + 2p)po + 48p"(32p+ + llp)p^ 



+60p(13p+ + 12p)pj^ + 9(9p+ + 35p)p^ 



(5.34) 
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Figure 5.14: Real and imaginary parts of the fundamental quasi-normal frequency for metric per- 
turbations (M = J = 0) \K\ = (blue), |i^| = 1 (red), \K\ = 2 (yellow) of squashed Kaluza-Klein 
black holes. Higher values of K correspond to higher oscillation frequency and slower damping. 



The fundamental quasi-normal modes calculated with 6th order WKB method (see sec. 3.3) 
are presented on the figure 5.14 The real part of uj decreases for higher values of po- The absolute 
value of imaginary part has maximum for K 0. For not large values of po the quasi-normal modes 
of such Kaluza-Klein black hole perturbations are longer lived. Their behavior for small values of 
Po in time-domain looks like massive field tails (see chapter [6]), which occur at much earlier time. 

It is important to note, that if we know the quasi-normal frequency we can find the size of the 
extra dimension in the considered black hole model, so that quasi-normal modes give a kind of 
opportunity to "look into" an extra dimension at low energies. In detail, when a dominant quasi- 
normal mode is measured, one can compare it with the numerically found one and find out which 
is the value of po and the radius of the event horizon that corresponds to the observed quasinormal 
mode. In this way we can determine the parameters of the black hole and the size of the extra 
dimension, assuming that there exist no other Kaluza-Klein black holes with similar features. Also, 
we do not observe any growing mode, what supports stability of such black holes. 
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Chapter 6 

Massive fields around black holes 



6.1 Evolution of massive fields. Quasi-resonances 



In the previous chapters we have considered massless fields propagating in the background of various 
black holes. Let us go further and study evolution of massive fields. 

One should note, that a scalar field with the mass term can be interpreted as a self-interacting 
massless scalar field within regime of small perturbations |101] . A massless scalar field, when 
considered in models with extra dimensions of Randall-Sundrum type [771 EH] , gains a large effective 
mass due to the Kaluza-Klein momentum. Also, the effective mass is acquired by a massless scalar 
field in the vicinity of the magnetised black holes |102j . 

Let us start from the consideration of a test massive scalar field in the Z)-dimensional Schwarz- 



schild background, given by the metric (5.1) with 



fir) 



r 



D-3 



The Klein-Gordon equation. 



d 



..ab 



-g dx°- 



9 



-g- 



'fl^{x) 



(6.1) 



which governs the massive scalar field evolution in the curved background, can be reduced to the 

3 potential 

D-3) {D-A){D-2) 



wave-like equation (2.10) with the effective potential 



V{r) = f{r) fi' + 



4r2 



-fir) + ^/'(r) 



2r 



(6.2) 



where fi is the field mass and / = 0, 1, 2 . . . is the multipole number. 



behavior of the eigenfunction (2.10), so that 



One can see, that the potential (6.2) does not vanish at the spatial infinity. It changes the 



+00 : oc exp(ixn 



(6.3) 



where x = a/^^^^ — /i^ (see sec. 3.4.7) 



Repeating the calculations of the section 2.3.2, instead of (2.27) we find [12 



Reix)\QUn = oo)|2 + Reiu)\Q^in = -oo)|2 + 2Re(u;)Im(cu) / \QUn)\^dn = 0. 



(6.4) 
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Since Re(x) has the same sign as Re(u;), the unstable modes do not oscillate (see sec. 6.4). 

For some values of the black hole mass M and the scalar field mass /x purely real frequencies 
were observed. Such oscillations have infinitely long lifetime and were called quasi-resonances [FF] . 
Indeed, for massive fields quasi-normal frequencies are not required to be complex. Let us suppose 
that for some parameters we have Im(u;) = 0, so that Re(x) = 0, implying 

Q^(r, = -oo) = 0. (6.5) 

It means, that there is no wave at the event horizon and, because Re(x) = 0, no energy transmission 
to the spatial infinity. That is why such oscillations do not decay. The situation is similar to the 
standing waves on a fixed string. The requirement for x to be imaginary, bounds the purely real 
frequencies by the field mass, what was pointed out in [T2] . 



^QRM < At- (6.6) 
The phenomenon of quasi-resonance exists since the effective potential is not zero at the spatial 



infinity. If the potential vanishes at — > ±oo the condition luiiu) = does not satisfy (2.27). 
Thus, for instance, there is no quasi-resonances in the spectrum of the massive scalar field in the 
Schwarzschild-de Sitter background. 



6.2 Quasi- normal spectrum of the massive scalar field 
around Schwarzschild black holes 

The quasi-normal modes of the massive scalar field around 4-dimensional Schwarzschild, Kerr |lU3j 
and Reissner- Nordstrom black holes [TD] and of the massive Dirac field in the 4-dimensional Schwar- 
zschild background |104j were calculated within the 3rd order WKB approach. Yet, the WKB for- 
mula, used for the calculations, is valid only for small field mass. For large mass of the field the 



WKB approach needs modifications (see sec. 3.3), that were not taken into account in the earlier 
research [701 [TU31 [TUl] . 

The quasi-normal spectrum of the massive scalar field for Reissner-Nordstrom black holes were 
studied for the first time by the more accurate Frobenius method by Ohashi and Sakagami [57]. 
They have found the quasi-resonances and supposed that they appear in a limiting situation and 
quasi-normal modes can disappear when the field mass exceeds a certain value. We found [I2] that 
after reaching a quasi-resonance only one quasi-normal frequency disappears. The higher overtones 
remain in the spectrum for any finite field mass. 

The quasi-normal spectrum of massive scalar fields in D-dimensional Schwarzschild background 
was studied for the first time in HHI within Frobenius method. 



On the figure QA_ one can see, that increasing of the field mass gives rise to decreasing of the 
imaginary part of the quasi-normal mode until reaching the vanishing damping rate. When some 
threshold values of /ir_|_ are exceeded, the particular quasi-normal modes "disappear". The larger 
field mass is, the more first overtones share this destiny. The disappearing of a finite number of 
modes implies that for large field mass some higher overtone becomes a fundamental frequency, 
having the longest lifetime among the rest of the modes. Nevertheless, the quasi-normal spectrum 
remains infinite for any finite field mass. 

The behavior of the fundamental mode depends qualitatively on the number D. For = 4, 5, 
the fundamental mode behaves like higher overtones: decreases its imaginary part, reaches quasi- 
resonance for some value of /i, and disappears. For > 6 we observe qualitatively different picture 
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Figure 6.1: Three higher quasi-normal modes (Z = 0) of the Schwarzschild black hole {D = 6) for 
the massive scalar field of various /x. The frequency for = has the largest imaginary part. The 
points were plotted with the step of Ayur+ = 1/2. Solid lines mark the same overtone number. 



see fig. 6.2). As fi grows the imaginary part of the fundamental mode tends to zero asymptotically 



while the real part approaches /i (see fig. 6.3) 



This property of the fundamental frequency behavior for large D leads to another remarkable 
fact. Since for high field mass the imaginary part of the fundamental mode tends to zero only 
asymptotically (see fig. 6.2), and the imaginary part of higher overtones reaches zero for some 

for which the imaginary parts of 



6.3), there are some values of fir^ 



(see fig. 6. 
finite value of /ir+ (see fig 
two overtones are the same. After one of these values is reached the overtones can be distinguished 
only by their real part and the quasi-normal ringing has two dominant frequencies in its spectrum. 
Thus one could observe the superposition of that two frequencies at late times of the quasi-normal 
ringing. 

The dependance on the field mass is qualitatively the same for all multipole numbers (see fig. 
6.3). Since the real part of the quasi-normal frequency of higher / is larger, the quasi-resonances 



are reached for larger n in order to satisfy the constraint (6.6). 
It was shown both analytically and numerically 



that asymptotically high overtones do 
not depend on the field mass and satisfy the same formula as for the massless case 



oon = (± ln(3) - 2m (n+^ 



THawking ( ± ln(3) - 2ni 



n 



It is important to note, that despite in most cases quasi-normal modes of the fields of different 
spin show qualitatively the same behavior, when the field is massive, the field spin becomes crucial. 
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Figure 6.2: Fundamental quasi-normal frequencies for the D-dimensional Schwarzscliild black hole 
for the massive scalar field (/ = 0) of various fi. The frequency for n = has the largest imaginary 
part. The points were plotted with the step of A/ir+ = 1/10. Solid lines mark the same number of 
D. 



The quasi-normal modes of the massive vector field for 4-dimensional black holes were studied in 
|105j . It was found that the fundamental mode shows correlation with the field mass totally different 
from all the remaining higher overtones. For massive scalar and massive vector fields the behavior 
of the fundamental frequency is qualitatively similar: the real part grows with fi, while the absolute 
value of the imaginary part decreases until reaching zero (quasi-resonance) and then disappears. 
The higher overtones have their real part decreasing to tiny values, and, the absolute value of the 
imaginary part is growing with /i, leading to existence of almost pure imaginary modes which just 
damp without oscillations. 



6.3 Stability and quasi-normal modes of the massive scalar 
field around Kerr black holes 

Analysis of the quasi-normal spectrum allows also to prove stability of massive fields in the Kerr 
background numerically [H]. The stability was proved analytically only for high field mass |106] : 

^ \m\a f 2M ^ 
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Figure 6.3: Fundamental quasi-normal frequencies of the Schwarzschild black hole {D = 6) for the 
massive scalar field as function of the field mass fi for I = 0, 1, 2 (represented as dots, rhombuses 
and triangles respectively). The dashed line corresponds to Re(u;) = /x. 



Strictly speaking, because the azimuthal number appears in the righthand side, for any finite field 
mass one could find such kind of perturbations, that is not proven to be stable. In practice, the 
numerical proof is reduced to checking that there is no growing mode in the quasi-normal spectrum 
for small field mass and high enough value of m [14]. 

The Kerr metric is described by (3.53 ), if we put Q = 0, A = 0. The Klein- Gordon equation for 
the massive scalar field (6.1) in this background, after substituting the ansatz 



iuit im4> 



= e e 

allows to separate variables. The radial part satisfies 



Rir)Sie), 



d 

dr 



A. 



dR{r) 



dr 



A. 



-A 



2 2 



R{r) = 0, 



where A is the separation constant 

A, 

The angular part 



2 I 2 



2Mr, 



K = cj(r2 + a') 



am. 



1 d 
sin 6 dd 



ds{e)\ 
de J 



m 



(6.8) 



(6.9) 



sin" 9 



- a'uj' sin^ 6 - cos^ 6 + 2amuj + A S{e) = (6.10) 
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Figure 6.4: High overtones of the Schwarzschild black hole for the massive scalar field {D = 4, 
/ = 0) for /ir+ = 0.6 (dots) and /ir_|_ = 6.0 (boxes). 



is equivalent to (5.13) up to the redefinition of variables and can be solved as described in the 
section 13.4.51 



The appropriate Frobenius series for (6.8) are given by 



R{r) = exp(ixr) 



2 \ ^f^-i 



rr I 



rr+ — 



A:=0 



rr I 



rr+ — 



(6.11) 



where x = \/ 



The coefficients of (6.11) satisfy the three-term recurrence relation (3.41), which allows to 



find the spectrum numerically and prove, that it does not contain unstable modes jl4j . 

We can see that despite the lifetime of the quasi-normal oscillation increases monotonously as 



the azimuthal number m grows, its imaginary part remains bounded (see fig. 6.5) and, thereby, 



does not show any tendency to instability. Therefore, we conclude that the formula (6.7) implies 
stability for the large field mass. 

The dependance of the quasi-normal frequency on the field mass fi is qualitatively the same as 
for the Schwarzschild case: larger fi leads to considerable decreasing of the damping rate and finally 



to the quasi-resonances (see sec. 6.2). 

The dependance on the rotation parameter a is more complicated. For m = the imaginary part 
of quasi-normal frequencies decreases its absolute value as a grows, while the real part increases, 
reaches maximum and then decreases. 
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Figure 6.5: The fundamental quasi-normal frequency of Kerr black holes for the massive scalar field 
as a function of azimuthal number m for I = Q, a = 0.15r+, /ir+ = 0.2. 



6.4 Quasi-normal modes of black strings and the Gregory- 
Laflamme instability 

Unlike four dimensional Einstein gravity, which allows existence of black holes, higher dimensional 
theories, such as the brane-world scenarios and string theory, allow existence of a number of "black" 
objects: higher dimensional black holes, black strings and branes, black rings and saturns and others. 
In higher than four dimensions we lack the uniqueness theorem, so that stability may be the criteria 
which will select physical solutions among this variety of solutions. 

According to the brane-world scenarios, if the matter localised on the brane undergoes grav- 
itational collapse, a black hole with the horizon extended to the transverse extra direction will 
form. This object looks like a black hole on the brane, but is, in fact, a black string in the full 
D-dimensional theory. Such black strings suffer from the so-called Gregory-Lafiamme instability 
[1U7[ I1U8] . which is the long- wavelength gravitational instability of the scalar type of the metric 
perturbations. The threshold values of the wave vector k at which the instability appears were 
found in [Tn9] . 

The evolution of the spherically symmetric linear perturbations of D-dimensional black strings in 
time and frequency domains was studied in p6] . For the first time the quasi- normal modes and time- 
domain profiles were studied in the stable sector. Also the appearance of the Gregory-Lafiamme 
instability was shown in the time domain. 
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Let us note, that the equations for gravitational perturbations contain the Kaluza-Klein mo- 
mentum. The contribution of this momentum in the corresponding effective potential looks like the 
mass term. 

For the static black string in D > 5 space-time dimensions, the background metric can be written 

as 

ds^ = gf.^dx^dx" = f{r)df - — - r'^dQl_.^ - dz^, (6.12) 

where ^ 

/(r) = l-(^)", n = D~A. 

and (if2|,_3 is the line element on a unit {D — 3)-sphere. 

The ^-direction is periodically identified by the relation z = z + 2'kR. Let us study the {D — 3)- 
spherically symmetric perturbations, which we can write in the following form 

771 

59,,u = e''''a^^{t, r), k = —, m E Z. 

ri 

The perturbed vacuum Einstein equations have the form 

5R^, = Q. (6.13) 
The perturbations can be reduced to the form, where the only non-vanishing components of a^y are 

dll Clj.^ /ij., O'zz ^zi (^tr ^vj (^zr ikfly. 

The linearised Einstein equations give a set of coupled equations determining the four radial 
profiles above. However, we may eliminate hy, and ht from these equations in order to produce 
a single second order equation for hz'. 

hz = f{rfK+p{r)h', + q{r)hz, (6.14) 



where 



p{r) = ^^^^ 1 + ^ ^ 



q{r) = -k'^fir) 



2A;V2 - n(n + 3) (r+/r)" 
2p^2 _^ ^(-^ _^ l)(r+/r)"' 



Defining 



^-(n-l)/2 

hz(t,r) = — TTrr; ; TT -— \I^(t,r), 

^' ' 2Pr2 + n(n + l)(r+/r)" ^' ^' 



we can reduce the equation (6.14) to the wave-like equation 



dv 

where dr^, = — — — is the tortoise coordinate. Here, the effective potential V(r) is given by 
f{r) 

V(r) = ^(") ^(") 

4r2 (2A;2r2 + n(n + l)(r+/r)")^ ' 
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Figure 6.6: Time-domain profiles of black string perturbations for kr+ = 2.5: n = 2 (red, top), 
n = 3 (orange), n = 4 (green), n = 5 (blue, bottom). Late-time decay of perturbations for n > 3 is 



where 



U{r) 



16^6^6 ^ 4^V(^^ ^ 5)(3/(r) - 2n + 3n/(r)) - 
-4:k'^r^n{n + 1) (n(n + 5) + /(r)(2n^ + 7n + 9)) 

-n^(n + 1)3 (/(r) - 2n + n/(r)) (^^j 



r 



2n 



One can see that k plays the role of the effective mass. At asymptotically late time we observe 



power-law damped tails, which have the oscillation frequency equal to k (see fig. 6.6), resembling 
asymptotical behavior of massive fields near Schwarzschild black holes (see sec. 6.5). The behavior 
of the first overtone is qualitatively similar to that of the fundamental mode for massive fields of 



higher- dimensional Schwarzschild black holes (see sec. 6.2). 

• For D = 5 {n = 1), as k grows, the imaginary part of the first overtone quickly decreases and 
vanishes for some threshold value of k, while its real part stays smaller than the threshold value. 
After the threshold value of k is reached, the first overtone "disappears" . 

• For D > 6 (n > 2), the imaginary part of the first overtone becomes small for large k, while 
the real part asymptotically approaches k. 

Even though the first overtone of the spherically symmetric black strings behaves similarly to the 
fundamental mode of massive fields near higher-dimensional Schwarzschild black holes, the other 
modes have completely different behavior. The fundamental mode of black string perturbations 
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Figure 6.7: Time-domain profiles of black string perturbations for n = 1 kr^ = 0.84 (magenta, 
top), kr+ = 0.87 (red), /cr+ = 0.88 (orange), fcr+ = 0.9 (green), fcr+ = 1.1 (blue, bottom). We can 
see two concurrent modes: for large k the oscillating one dominates, near the critical value of k the 
dominant mode does not oscillate, for unstable values of k the dominant mode grows. The plot is 
logarithmic, so that straight lines correspond to an exponential growth or decay. 



is purely imaginary. It grows for small values of k, leading to instability of the black string (see 



fig. 6.7). At moderately large values of k, sufficiently far from instability, the profile has the same 



form as that for massive fields, yet, when approaching the instability point, the real oscillation 
frequency and the decay rate decrease considerably. After crossing the instability point we observe, 
that starting from some tiny values, Im(a;) > are slowly increasing (while Re(a;) = 0). Therefore 
we can conclude, that there is some static solution u = of the wave equation, which shows itself 
exactly in the threshold point of instability. 



6.5 Late-time tails of massive fields 

The late-time behavior of black hole perturbations was studied for the first time by R. Price, 
who showed that perturbations of the massless scalar and gravitational fields decay as oc t"^^'"'"^) 
at asymptotically late time |110[ llllj . In [112j Bicak found that the scalar massless field in the 
Reissner- Nordstrom background decays as oc ^-(^'+2) \Q\ < M and as oc for the extremal 

black hole charge \Q\ = M. For Schwarzschild-de Sitter and Reissner-Nordstrom-de Sitter black 
holes, instead of power-law tails, the exponential tails were found |36j . 

In higher dimensional space-times the late-time behavior depends also on the number of extra 



74 



dimensions. It was found that massless scalar and vector fields, and gravitational perturbations of 
higher- dimensional Schwarzschild black holes have the decay law oc j(~(2'+£'-2) jj yn^^ 

oc t-(2'+30-8) ^Qj, even D > A |114j . The same late-time behavior was observed for odd- dimensional 
Gauss-Bonnet black holes [89] . The late-time tails of the brane localised Standard Model fields were 
studied in P5|- Their late-time decay law is oc t-C^'+^^-i)^ 

One should note that the late-time decay of perturbations within the full non-linear gravity does 
not agree with the linearised theory in dimensions higher than four. If we consider a massless scalar 
field and take into account the back reaction of the field upon the metric, the late-time decay rate 
becomes smaller, coinciding with the linearised theory prediction only in four dimensions |115] . 

The late-time behavior of massive fields is qualitatively different from massless ones: at late 
time the decay profile is oscillatory inverse power law tail. Also, the field mass implies different 
behavior at the intermediate late time (1 -C t/M < {^M)'^) and at the asymptotically late time 
(t/M > [iiM)-^). 

For a massive scalar field with mass /i in the background of the Schwarzschild black hole, the 
perturbations decay as oc t~('+3/2) sin(/it) at intermediate late time [10 Ij and as oc 

t-(5/6) sin(^i) at 

asymptotically late time [1161 1117j . The same behavior at asymptotically late time was found also 
for the massive scalar field perturbations of the dilaton black hole [118j and the Kerr black hole 
[119] and for the massive Dirac field in the Schwarzschild background [72] . 

For higher dimensional Schwarzschild black holes the intermediate late-time behavior of the 
massive scalar field is found to be oc sin(/it) [120] . 

The late-time behavior of the massive vector field was studied both numerically and analytically 
for the first time in [T5]. At intermediate late time the decay law depends on the polarisation, being 
either oc t^^'+^Z^) sin(/xt), or oc t^C+^Z^) sin(/it), or oc sin(/it). At the asymptotical late time 

the behavior is the same, as for other massive fields, oc sin(/it). Therefore, we conclude that 

the asymptotically late-time decay law does not depend on the spin of the field. 

The late-time decay law for massless perturbations, the intermediate and asymptotically late- 
time behavior of massive field perturbations for various spins in asymptotically fiat backgrounds 
are presented in the following table. 



perturbations 


massless 


massive (intermediate) 


massive (asymptot.) 


s = (4D Schw.) 
s = i \aD Schw.) 
s = 2 (4D Schw.) 
s = (4D Kerr) 
s = (4D R-N) 


^-{2«+3) 
^-(2«+3) 
^-{2«+3) 
^-(2«+3) 
■f--{2l+2) 


t-('+3/2) gin(^^) 
{r('+l/2),t-('+3/2),t-('+5/2)|sin(^i) 

t-{'+3/2) sin(^t) 


t-(5/6) sin(/it) 
t-(5/6) sinlut) 

t-(5/6) sin(^t) 


s = (odd D) 
s = (even D) 
on a 4-brane 


^-{2l+D-2) 
^-(2l+3D-8) 
^-{2l+D-l) 


tHl+iD-l)/2) 

t-('+(^-i)/2)sin(/xi) 
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Chapter 7 

Quasi-normal modes of black holes, 
whose metrics are unknown analytically 



7.1 Numerical methods 



Black holes are compact objects. Therefore, it is natural to expect that their quasi-normal ringing, 
being the property of the black hole, at least in the dominant order, does not depend on what 
happens at large distance from the black hole. 

It turns out that the dominating frequencies depend mostly on the black hole solution behavior 
in some region near the event horizon. Thus one can find them in frequency domain, even if the 
behavior of the solution at large distance is not known. 

Two qualitatively different examples were considered in this context: the scalar hairy black 
hole in the anti-de Sitter background and the Einstein- Aether black hole in the asymptotically flat 
background. In both types of the background the behavior of the solution at large distance is not 
important [28] : 

• For an asymptotically anti-de Sitter background we usually require Dirichlet boundary condi- 
tions at spatial infinity. The most significant part of the metric perturbations stays, thereby, near 
the black hole. That is why the solution behavior at this region causes dominant influence on the 
quasi-normal spectrum [13j. 

• For the asymptotically flat case the searching of the quasi-normal modes can be reduced to 
the scattering problem. Therefore, the quasi-normal frequencies are determined mainly by the form 
of the effective potential near its peak [121 US] . 

In order to calculate quasi- normal frequencies for the asymptotically anti-de Sitter background 



one can use the Horowitz-Hubeny method, described in the section 3.5 It is clear that if s{z), t{z) 



and u{z) in (3.71) are series which converge at the spatial infinity {z = 1) quickly enough, then yn 
is still possible to calculate because it depends insignificantly on higher terms of the series. In order 
to find series expansion for s{z), t{z) and u{z), one can use equations which define the black hole 
solution. We can always do this because s{z), t{z) and u{z) can be explicitly expressed in terms of 
the metric coefficients and their derivatives. 

Thus, the quasi-normal modes could be found without solving the equations for the black hole. 
It is enough to find series expansions for the metric coefficients near the horizon. It is important 
to note that one can control the precision of the eigenfrequencies u by requiring the convergence of 
the found result with respect to increasing of the number of expansion terms of all the series [13j . 

For asymptotically fiat black hole solutions it is convenient to use the WKB method (see sec. 



3.3), where the asymptotic solutions of the wave equation near the event horizon and near spatial 
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Figure 7.1: The effective potential for electromagnetic perturbations near the Schwarzschild black 
hole (r_|_ = 1, Z = 2) and the same potential interpolated numerically near its maximum. 



infinity are matched with Taylor expansion near the peak of the potential. We conclude therefore 
that the low-damping quasi-normal modes are determined mainly by the behavior of the effective 
potential near its peak. 

This statement was checked by considering the potential for the Schwarzschild black hole and 
also two other potentials, which lay closely to the Schwarzschild potential near its maximum, but 
have very different behavior far from the black hole. These two potentials are chosen in the following 
way. We choose some points near the maximum of the effective potential. The first potential is an 



interpolation of these points by cubic sphnes (see fig. 7.1). The second potential is a fit for the 
points by a ratio of polynomial functions. 

Since the WKB formula contains the value of the effective potential and its derivatives at the 
potential peak, we find that the results obtained with the help of all three potentials lay very close 
if the interpolation and the fit were made with the appropriate precision |16|. Despite the higher 
derivatives of our interpolation potential are not defined, we are able to evaluate them step by 
step by interpolating in the same way the first and all the consequent derivatives of the potential. 
It turns out that the interpolation potential is very sensitive to numerical errors. Therefore, to 
calculate the quasi-normal modes with the appropriate precision, one must calculate the values of 
the potential with very high accuracy. In fact, for the practical purposes one can use fitting of the 
potential which does not accumulate the numerical error. 

To test the accuracy of this approach one can use the fact that the sixth order WKB formula gives 



a smaller relative error than the third order one (see sec. 3.3). Since the higher WKB order depends 



on the higher derivatives of the effective potential, that are more sensitive to the interpolation or 



78 



fitting accuracy, the higher order WKB formula should give some random values, if the accuracy is 
not enough. 

The independence of the quasi-normal modes on the behavior of the effective potential at large 
distance, being inspired by the WKB formula, is not related with the WKB method. However, since 
the WKB formula depends only on the potential near its peak, we unable to prove the statement 
within the described approach. That is why, it is important to check the results with a different 
method. In order to do this, we have found in time domain the quasi-normal modes of gravitational 
perturbations of black holes in the Einstein-Aether theory. The calculated frequencies show an 
excellent agreement with the WKB results |19j . 



7.2 Quasi- normal modes of the scalar hairy black hole 

Since the paper of J. Bekenstein [121j . it is well-known that black hole can not have scalar hair 
within minimal coupling. However, as have been found in |122l 1123] , there is a possibility of dressing 
a four-dimensional black hole in anti-de Sitter space-time with a non-minimally coupled classical 
scalar field. In |124j the results were extended for higher dimensional configurations. Despite we live 
in the de Sitter universe, black holes in the anti-de Sitter background attracted considerable interest 
due to the AdS/CFT correspondence In this context, the quasi-normal ringing of massive non- 
minimally coupled scalar field to the black hole in the anti-de Sitter background was studied for the 
first time in [13]. 

We consider the spherically symmetric solution 

ds^ = N{r)e^^^'^Me - N{r)-'dr^ - [dO^ + sin 9^ V) (7-1) 

of the action, which describes a self-interacting scalar field with non-minimal coupling to gravity: 

'I/O OAA /"V" 



S = d X \/—g 



{R-2K)--{V4>Y--iR<P' 



(7.2) 



_2 ' ^ 2 ' 2" " 2 

where R is the Ricci scalar, A is the cosmological constant, ^ is the coupling constant, and 

(V0)' = Va0V>. 

In order to simplify the equations of motion we use the following conformal transformation 

gal, = (1 - e0')^7af,. (7.3) 

After the transformation, the action takes the form 



S = j d'x^g(^- {R-2k)-\{V^f-Um^ 



(7.4) 



where we define 



, .O^C^mW jj(<,)^'±±3ll^ytL± (7 5) 



The metric takes the following form 

df = N{f)e^~^^'Ut' - N{f)-'df^ - f2 {de^ + sin 6^ dip"") , (7.6) 
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where 



N = iV(l-e0'-er00')'(l-e0')"', 
Varying the action, one can find the following equations of motion [123] 

2 



d{fN) 
df 



. A 2 2 / ^ f 



d6 


r / d^ 


df 


2 \ df 









= N 


dr'^ 



2 Kdf) 



rfA^ 2N\ d^ dU 
df df f I df d^ 



(7.7a) 

(7.7b) 
(7.7c) 



In order to obtain the series expansion for the metric coefficients (7.6), we have to impose 
some boundary conditions, which are associated with the black hole parameters. We are using the 
following parameters. The event horizon radius f_|_ can be chosen arbitrary in order to fix length 
scale. This choice defines the boundary condition for N: N{f+) = 0. In order to measure all 
dimensional values in units of f+, we choose f_|_ = 1. The value S{f^) is chosen arbitrary in order to 
fix time scale. In order to introduce the boundary conditions in the same point we impose 5(r+) = 0, 
what differs from Winstanley's S{oo) = 0. The last parameter (f){f+) = Q can be associated with 
the scalar charge of the black hole. 

Spherically symmetric perturbation equation for 



(5$ 



takes the wave-like form (2.21), where the perturbation potential is given by |122j : 



V{f) 



Ne 



25 



(f/ + A)r-)|l-r- (f)' 



N + 2f 



dU d<^ 
d^ df 



f2 



d^U 

d$2 



(7i 



The tortoise coordinate is given by 



dr* 
df 



It is similar to the ordinary anti-de Sitter case: 



(7.9) 



r = r_(_ 
f = oo 



oo, 
0. 



This perturbation corresponds to infinitesimal changing of the black hole mass and is similar to 
the zero-multipole scalar perturbations near a hairless black hole. The crucial difference is that the 
perturbation of a scalar field due to the hair gives a first-order correction to the metric. 

In order to calculate the quasi-normal modes we find the series expansions near the black hole 
event horizon for $, 6 and N using the equations (7.7). Since N grows at the spatial infinity, it is 
convenient to find the series for the function g = N /f"^ instead of N. After the series for the metric 
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Figure 7.2: Dependance of the real and imaginary parts of the quasi-normal frequency 
on Q for the scalar hairy black hole. The imaginary part quickly decreases as the "scalar 
charge" Q increases. The real part reaches its maximum and then quickly falls down to zero. It 
vanishes for some non-critical charge Qq < For Q > Qo the real part remains zero and the 

frequency is purely imaginary (within numerical precision). 



functions are known one can find the series for the effective potential (7.8) and, finally, for s{z 



t{z) and u{z) in (3.71). Then we use the technique, described in the previous section. 

It was found [13] that imaginary part of the quasi-normal modes changes monotonously as a 



function of A, Q and ^. Scalar field mass fi cause imaginary part to decrease its absolute value but 
not to vanish. It reaches a minimum for some particular mass and then increases very rapidly (see 



fig. 7.5). Thus no infinitely long living oscillations appear. 



Real part has more complicated behavior and depends significantly on the parameters Q and A 



see figures 7.2 and 7.3). Other parameters, ^ and fx, change the real part within the comparable 



small bounds (see figures 7.4 and 7.5). For relatively large values of ^, A and Q purely imaginary 
frequencies appear. 



From the figure |7.6| we see that higher overtones tend to equidistant spacing similar to hairless 
black hole spectrum behavior in the anti-de Sitter background |125j . 



7.3 Perturbations and quasi-normal modes of black holes 
in the Einstein- Aether theory 

General relativity is based on the local Lorentz invariance. Yet, there appeared a lot of attempts to 
go beyond the local Lorentz symmetry. Aether can be considered as locally preferred state of rest at 
each point of space-time due-to some unknown physics. Einstein-Aether theory is general relativity 
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Figure 7.3: Dependance of the real and imaginary parts of the quasi-normal frequency 
on A for the scalar hairy black hole. The real part increases for relatively small "scalar charges" 
Q, but for large Q = 2.2 it also reaches maximum and then falls down and vanishes for A ^ —0.292. 
For larger values of A the real part remains zero and purely imaginary frequencies exist. The 
imaginary part decreases almost linearly as A grows. 



coupled to a dynamical time-like vector field, which is called "aether" . This theory is what comes 
instead of usual general relativity when local Lorentz symmetry is broken. Namely, the vector field 
breaks local boost invariance, while rotational symmetry in the preferred frame is preserved (see [1| 
for review). 

It is important that the significant difference between Einstein and Einstein-Aether theories 
should be seen in the regime of strong field, for instance in observing of the characteristic quasi- 
normal spectrum of black holes. Thus, existence of aether could be tested in the forthcoming 
experiments with new generation of gravitational antennas. In this context, the quasi-normal spec- 
trum for test scalar and electromagnetic fields [16] and for the gravitational perturbations [19] was 
studied. 

The solution for a black hole metric in the Einstein-Aether theory was found numerically in 
|126[ 1127] . Since the analytical solution is not known, in order to calculate the quasi-normal modes. 



one can apply the fitting of the effective potential as described in the section 7.1 



The Lagrangian of the full Einstein- Aether theory that forms the most general diffeomorphism 
invariant action of the space-time metric Qab and the aether field u"", involving no more than two 
derivatives, is given by 



-R - K''\^Vau"'Vbu'' - XigabU^u' - 1), 



(7.10) 
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Figure 7.4: Dependance of the real and imaginary parts of the quasi-normal frequency 

on ^ < ^ for the scalar hairy black hole. The imaginary part decreases as ^ grows. The 
quickness of such decreasing depends on A and Q. For some A and Q (as presented on the bottom 
graphic) the real part can fall down to zero and purely imaginary frequencies appear for larger ^. 



where R is the Ricci scalar, A is a Lagrange multiplier, which provides the unit time-like constraint. 



imni 



where Cj are dimensionless constants. 

Spherical symmetry allows to fix C4 = 0. Following |126l 1127] . we shall consider the so-called 
non-reduced Einstein- Aether theory, for which C3 = 0, and we can use the field redefinition that 
fixes the coefficient C2: 



C2 



2 - 4ci + 3cf 



,2 ' 



so that Ci is the free parameter. 

The metric for a spherically symmetric static black hole is 

ds"" = N{r)de - ^^dr^ - r^d^f , (7.11) 

where the functions N[r) and B[r) are given by numerical integration near the black hole event 
horizon. 



The perturbation equations can be reduced to the wave-like form (2.10) with the effective po 
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Figure 7.5: Dependance of the real and imaginary parts of the quasi-normal frequency 
on /i for the scalar hairy black hole. We consider only < /i^ < — 4^A. The real part changes 
within comparatively small range. The imaginary part reaches its maximum remaining negative. 
For this value of /i the oscillations have the longest lifetime, but they still damp and we do not 
observe quasi-resonances. For larger scalar field masses the imaginary part of frequency decreases 
and the damping rate is higher. 



tentials 



+ , 1 d N{r) 
r dr^ B{r) 



N{ry^^^ + ~^^, / = 0,1,2, 

([+2)([^ ^ 2N\r) _ 1 d N{r)^ 
J.2 B^(r)r'^ r dr^ B(r) 



I = 2,3,4 



(7.12a) 
(7.12b) 
(7.12c) 



for scalar. Maxwell and axial gravitational perturbations respectively. 
The tortoise coordinate is defined as 



dn 



B{r) 
iV(^ 



-dr. 



Note, that since the background value of aether coupling is small in comparison with the back- 
ground characteristics of large black hole, the Schwarzschild metric is slightly corrected by the 
aether. That is why, when considering perturbations of the metric we neglect small perturbations 
of aether, keeping only linear perturbations of the Ricci tensor. 

Since N{r) and B{r) are known numerically, we can use the fitting technique [16]. The fit 
functions for N{r) and B{r) were chosen as fractions of two polynomials, which are characterised 



84 



Re OJ 
2 - 



1.5- 



1 - 



. 5 



^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ -Imo) 

0.1 0.2 0.3 0.4 0.5 0.6 

Figure 7.6: First five quasi-normal frequencies for the scalar hairy black hole A = —0.1, 
Q = 1, C, = 0.1, /i = 0. The high overtones approach equidistant spacing. 



by the number of terms in their numerators and denominators. There is some optimal number for 
which the convergence of WKB series is best. Practically, in order to find optimal number of terms 
for N{r) we search for minimal difference between third and sixth order WKB values, first for the 
case B{r) = 1. When we have found the optimal fit for N{r), in a similar fashion, i.e., by looking for 
best WKB convergence, we are in position to find the optimal fit for B{r). Quick WKB convergence 
shows that higher derivatives of the metric coefiicients are calculated with the best accuracy. 

We find that as ci grows, the fundamental quasi-normal frequency increases the absolute value 



of its real and imaginary parts. This was checked also in time-domain [19]. On the figure 7.7 we see, 
that for higher ci the oscillation period and the lifetime of the perturbations decrease. We conclude 
that aether, if it exists, could be indirectly observed through detection of characteristic spectrum 
of black holes. 
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Figure 7.7: Evolution of axial gravitational perturbations (/ = 2) in time domain for non-reduced 
Einstein- Aether theory, ci = 0.1 (green line) and ci = 0.4 (red line), in comparison with the Schwar- 
zschild case Ci = (blue line). The higher ci is the quicker decay of the observed perturbations. 
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Chapter 8 

Perturbations of Schwarzschild black 
holes in laboratories 



8.1 Acoustic analogue of gravity 

In addition to the possibility to observe quasi-normal modes of black holes with the help of a new 
generation of gravitational antennas, there is a window for observation of the acoustic analogue of 
a black hole in laboratories. This is the well-known Unruh analogue of black holes |128[ 1129] . which 
are the apparent horizons appearing in a fluid with a space-dependent velocity, in the presence 
of sonic points. The wave, which passed through the sonic point, cannot propagate backward, 
mimicking thereby, the effect of the horizon at sonic points. 

The dynamic of sound waves in a fluid can be described by the propagation of a scalar field in 
some effective background. It turns out, that steady spherically symmetric flow cannot reproduce 
the Schwarzschild geometry because of the equation of continuity (see |130] for review). One could 
construct either an acoustic analog, which is conformal to the Schwarzschild black hole. 



ds^ oc ((l - y) dt^ - (l - y) ' dr^ " r^dO^ + sm9^d(l)^)^ 



or an exact analog to the 7-dimensional black hole, projected on the 4-brane (see section 5.3) 

-1 



ds^= ^1- i^'^Yj dt^ - (^1- (^"^y j dr^-r\de^ + sme^d(f'). (8.2) 

If one had a complete analogy with some known solution of the Einstein equations, say, the 
Schwarzschild solution, he could see, in the acoustic experiments, not only qualitative, but also, up 
to an experimental accuracy, exact numerical coincidence with a prototype characteristics. Namely, 
for quasi-normal modes, which are governed by the form of the wave equation, this numerical 
correspondence would mean that the effective potential of the perturbations of some hydrodynamic 
system coincides with an effective potential of the black hole. Fortunately, the consideration of the 
perturbations of a gas in de Laval nozzle |131j gives us such an opportunity of finding a system that 
obeys the same effective potential as a Schwarzschild black hole does [2T] . 



8.2 de Laval nozzle 

The canonical de Laval nozzle is a convergent-divergent tube, narrow in the middle. It allows 
to accelerate the gas until the sonic speed in its throat, reaching supersonic speeds after passing 
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the throat. The perturbations of the gas in de Laval nozzle can be considered as one-dimensional 
if the section does not change too quickly along the length of the nozzle. Let us show that the 
corresponding effective potential for perturbations in a canonical de Laval nozzle can be made equal 
to the potential for perturbations of Schwarzschild black holes, when choosing some specific form 
of the nozzle. 

We assume that a gas in the nozzle can be described by the equations of motion for the perfect 
fluid and that the flow is quasi-one-dimensional. The equations of continuity, the momentum and 
energy conservation read respectively 



dt{pA) + d^{pvA)=Q, 
dt{pvA) + d,[{pv^ +p)A] = Q^ 
dt{tA) + d,[{t + p)vA] = Q. 



(8.3a) 
(8.3b) 
(8.3c) 



Here p is the density, v is the fluid velocity, p is the pressure, A is the cross section of the nozzle. 



and 



-pv^ + - 
2 7 



P 



U) 



is the energy density. The heat capacity ratio for di-atomic molecules of air is 

7 = 1 + 2/72 = 7/5 = 1.4 (n = 5). 
We shall assume that the flow has no entropy discontinuity. Then the fluid is isentropic 

p (X p^ . 



Instead of (8.3b), we can use the Euler equation 

p{dt + vd^)v = 



-dxp. 



(8.5) 



For isentropic fluid, (8.6) is reduced to the Bernoulli equation 

dt<^ + ^{d,<i>y + h{p) = o, 

where h{p) = J p^^dp is the specific enthalpy and $ = J vdx is the velocity potential. 
According to |131] . the perturbation equations in such a nozzle can be reduced to 



•7) 



d^ 



dx 



*2 



+ K^-V{X*] 



K 



V{x*) 



1 


' 9 








_ 2 dx*^ 




\dx* ) 



8.9) 
.10) 



Here c^o = \/lRT/ p is the stagnation sound speed, which is used to measure x* in meters. The 
variable x* is an acoustic analogue of the tortoise coordinate which satisfies 



X [X 



+oo) = +00, 



namely. 



X 



X [X 



c.dx 



0) 



-oo, 



(8.11) 
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The function if^ represents small perturbations of the gas flow, 



pA 



fix) 



The sound speed is given by 



|f I dx 
c1 — f 2 



dp 7p 



(8.12) 

(8.13) 
(8.14) 

(8.15) 



" dp p 

Since (8.8) is invariant with respect to the re-scaling of g, we can fix the coefficients in ( |8.13[ ) 
arbitrarily: 

(8.16) 



9 2p(^-i)/2' 

Up to a coefficient, A can be found in terms of p as [21] 

A-=(l-p(-^))^/%. 

We find 

9 



P 



(l-7)/2 



P 



(1-7) 



Hence it follows that 



2(1 -p(7-l))l/2 2 (p(l-7) - 1)1/2 • 

pi-^ = 2/ (l ± ^/T^) . 



(8.17) 



(8.19) 



The sign in (8.19) should be chosen in order that p be a monotonous function with respect to the 



transverse coordinate. As we will show later, the function g for the Schwarzschild black hole can 
be chosen also monotonous in the R region, finite at the horizon and infinite at the spatial infinity. 
Therefore, we choose the minus sign. 



P 



1-7 



2^2 _ , g>l. 



(8.20) 



Substituting (8.19) into (8.17), we find the cross-section area as a function of g, 

1/(7-1) 



A 



V2 2g^l-y^ 



9' 



(8.21) 



1 



1-9- 



For the steady isentropic flow, (8.7) can be rewritten as 

,2 



V 



4t (""^ - 1) = ^ (i - v^^) - 1 



ci 7 - i 7 
Since f = at the event horizon, g must be finite there, and 

7 + 1 3 



9 



e.h. 



> 1. 



(8.22) 



(8.23) 



This requirement fixes both constants of integration. 



89 



8.3 de Laval nozzle for the Schwarzschild black hole 



Since g allows to find the cross section A and, thereby, the form of de Laval nozzle, if we find such 



g(x) that leads to the same expression for the nozzle potential (8.8), as the effective potential for 
the Schwarzschild black hole [21j . 



y(r) 



/(/ + 1) 



1 



ir+ = 1), 



(8.24) 



where / > s is the multipole integer number. The integer s = 0, 1, 2 describes the perturbations of 
fields of different spin: s = for the test scalar field, s = 1 for the Maxwell field and s = 2 for the 
gravitational perturbations of axial type. 



In order to find g that produces the required potential (8.24), we identify the "tortoise" coordi- 
nates of the black hole solution and of de Laval nozzle 



dr* = dx* 



p('-^y^dx 



2g^ (l 



1 — )dx 



1-r 



(8.25) 



The equation (8.25) relates the real coordinate of the nozzle x and the radial coordinate of the 



Schwarzschild solution r. It is convenient because we are able to find the equation for g{r) explicitly, 

j{r)r{r)g\r) + ]{rfg"{r) jirfg'irf 



2g{r) 



4^(r) 



V{r) 



(8.26) 



This implies that the form of de Laval nozzle is parameterised by r and x*[r) = r*{r). Note, that 
as we chose the radius of the event horizon to be unity, the nozzle coordinate x is measured in the 
units of the radius of the event horizon. 



The general solution of the equation (8.26) contains two arbitrary constants. They can be fixed 



in a unique way by the condition (8.23). Namely, the requirement that the solution must be finite 



at r = 1 fixes one of the constants. Then the other constant re-scales g{r), and must be fixed by 



its value at r = 1. Finally, the solution of (8.26), for arbitrary / and s, that satisfies (8.23), is given 
by the following formula. 



9[r) 



7 + 1 



2V2Vl~- 

7 + 1 
27V7^ 



-!)"+"(/ + n)! 



{n + s)\{n - s)!(/ - n)\ 



r(l + / + g)2Fi(g - /, g + / + 1, 1 + 2g, r) 

r(i + /-s)r(i + 2s) 



1.27) 



One can easily check that the above solution indeed satisfies the equation (8.23), for any fixed / and 
s. 



The equation (8.25) allows to find the dependance of the transversal nozzle coordinate x on the 
parameter r, 

+ 1 _ 4g^rf (l - v/l-^7(r)-2)) dr 



X 



1 /(r)(7-l)^2^?(r)2(l-yr3^M^) 
The integration constant is chosen so that x vanishes at the sonic point. 
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Figure 8.1: The form of de Laval nozzles and the effective potential in the nozzle coordinates. 



Now we are in position to find the required form of de Laval nozzle, i.e. to find its cross- 
section A{x). We just need to replace g{r) given in (8.27) and go over to the transverse nozzle 
coordinate x. The radius of the de Laval nozzle as a function of the transverse coordinate x is 



shown in figure |8.1[ Note that the canonical de Laval nozzle is diverging at the end of the flow 
trajectory, so that Ax=oo = oo. Indeed, the formula (8.27) implies divergence at least as ~ r^. 
Nevertheless, the diverging of the nozzle still let us keep the one- dimensional representation of the 
motion, because x is measured in units of black hole radius, i.e. one can "pull" the nozzle along 
the transverse coordinate x in order to make the area of the nozzle change as slowly as one wishes. 
Such a "pulling" simply means that we are getting the correspondence with a larger black hole. 

In a similar fashion, one can find the de Laval nozzle form for the gravitational perturbations of 
polar type [21]. The form of the nozzles for modeling polar and axial gravitational perturbations 
are very slightly different. 

It should be recalled also that the precise acoustic analogy is only established for a scalar field. 
To be able to reproduce the potential V{r) for fields of different spins, certainly, does not mean that 
one can reproduce all the characteristics of those fields in an acoustic model. 

The obtained acoustic analogue for the perturbations of the Schwarzschild black holes is not 
limited by quasi-normal mode problems only, but allow general investigation of propagation of 
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fields, including such processes as scattering and tunneling of waves and particles. 
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Chapter 9 
Summary 



Let us summarise the results of this dissertation. 

• We have performed the complete study of the influence of the cosmological constant on the 
quasi-normal spectrum of the Schwarzschild black hole for fields of various spins. We have found that 
the presence of the cosmological constant decreases the absolute value of the real and imaginary 
parts of quasi-normal frequencies. Also we have found analytically the large multipole limit for 



quasi-normal frequencies of Schwarzschild-de Sitter black holes (sec. 4.1). 

• We have studied the behavior of high overtones for the massless Dirac and massive scalar 
fields in the Schwarzschild background and for the test electromagnetic field and gravitational 
perturbations of the Schwarzschild-de Sitter black hole. We have shown numerically that the real 
part of the quasi-normal frequency asymptotically approaches zero for the Dirac and electromagnetic 
fields, but has an oscillatory behavior for the gravitational perturbations of the Schwarzschild-de 
Sitter black hole. These results were confirmed later analytically (sec. 4.2). The behavior of high 
overtones for the massive scalar field was studied both analytically and numerically. We have shown 
that the asymptotical behavior does not depend on the mass of the field, coinciding with the known 



analytical formula for the massless scalar field (sec. 6.2). 

• We have investigated the quasi-normal spectrum of the electrically charged scalar and Dirac 
fields in the background of the Kerr-Newman-de Sitter black hole. Special attention was paid to 
the influence of the electromagnetic interaction between the black hole and the test field upon the 



quasi-normal spectrum (sec. 4.3) 



• We have studied the quasi-normal spectrum of gravitational perturbations of the neutral and 
electrically charged black holes in higher dimensions with a positive cosmological constant. We have 
considered all kinds of the perturbations and prove that higher dimensional Schwarzschild-de Sitter 
black holes are stable. We have shown that the higher dimensional black hole is unstable if both the 
black hole charge and the cosmological constant are large enough. We have found the parametrical 
region of the instability (sec. 5.1). 

• We have considered gravitational perturbations of higher dimensional black holes in the 
Einstein-Gauss-Bonnet theory. We have found that the presence of the Gauss-Bonnet parameter de- 
creases the imaginary part of quasi-normal frequencies, causing the oscillations be longer-lived. Also 
we have studied the developing of instability of Gauss-Bonnet black holes in five and six dimensions 
(sec. 



5.2). 



• We have found quasi-normal modes and late-time tails of the rotating black holes and the 
non-rotating Gauss-Bonnet black holes for the Standard Model fields localised on a 4-brane (sec. 



5.3D . 

• We have studied quasi-normal modes of the Kaluza-Klein black holes with squashed hori- 
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zons. We have considered the test scalar field near rotating squashed Kaluza-Klein black holes and 
gravitational perturbations of non-rotating squashed Kaluza-Klein black holes (sec. 



5.4). 



• We have calculated the quasi-normal spectra of the massive scalar field in the backgrounds of 
the Tangherlini black hole (sec. 6.2), the Kerr black hole (sec. 6.3) and the scalar hairy asymptot- 



ically anti-de Sitter black hole (sec. 7.2). We have found that the massive scalar field is stable in 



these backgrounds and provided a comprehensive discussions about properties of their quasi-normal 
spectrum and quasi-resonances (sec. 6.1). Also we have studied the late-time tails for the massive 



vector field in the Schwarzschild background and found that the asjTiiptotically late-time decay law 
does not depend on the spin of the massive field (sec. 



6.5). 



• We have found quasi-normal modes and tails of the gravitational perturbations of black strings. 
Also we have considered the developing of the long-wavelength instability of a black string in time 
domain 



sec. 



6.4). 



• We have studied the infiuence of the local Lorentz symmetry breaking within the Einstein- 
Aether theory on the quasi-normal spectrum of the Schwarzschild black hole. We have found that 
the presence of the aether increases the absolute value of the real and imaginary parts of the 
fundamental quasi-normal frequency (sec. 7.3). 

• We have proposed a possibility of observation of the acoustic analogue of the Schwarzschild 
black hole in a de Laval nozzle. We have found the particular forms of the de Laval nozzles for 
which the effective potentials will coincide with the effective potentials for test fields or gravitational 
perturbations of the Schwarzschild black hole (chapter [s]). 

• We have developed two new numerical tools: 

1) the generalisation of the NoUert improvement of the Frobenius method for higher dimensional 
problems, which provides better convergence of the numerical procedure (sec. 



3.4). 



2) the method for the calculation of the quasi-normal frequencies of a black hole, which metric is 
not known analytically, but can be found as a numerical solution of a set of differential equations 
(sec. 



7.1). 



All the results reported here were obtained by us for the first time, and are, thereby, new. 
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